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ABSTRACT 

A  digital  simulation  language  for  the  interactive  definition  and 
solution  of  piece-wise  continuous  non-linear  ordinary  differential 
equations  using  the  Runge-Kutta  method  has  been  designed  and  imple- 
mented.    The  combination  of  interactive  graphics  approach  and  a 
special  differential  equation  description  language  make  the  analysis 
program  very  versatile  and  easy  to  use.      For  second  order  systems, 
a  grid  of  phase  trajectory  segments  over  the  user  specified  phase 
plane  is  used  as  a  background  for  the  solutions. 
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I.       INTRODUCTION 

The  goal  of  this  study  was  to  develop  an  interactive  analysis  tool 
for  the  study  of  piece  wise  continuous  non-linear  ordinary  differential 
equations.      It  is  recognized  that  there  are  available  in  most  computer 
centers  several  programs  for  solving  non-linear  differential  equations. 
However  they  are  either  very  special  purpose  or  else  require  that  the 
problem  and/or  control  description  be  inserted  into  a  subroutine  or 
part  of  the  analysis  program.     This  in  turn  must  be  compiled  by  the 
computer  before  execution  time.     If  the  user  wants  to  change  some 
feature  of  the  system,   he  must  change  the  applicable  program  cards 
and  recompile.     Thus  in  general  only  one  problem  type  per  computer 
run  is  possible. 

Also  there  exists  a  large  number  of  continuous  digital  simulation 
languages  [Clancy  and  Fineberg,    1965].     These  are  best  represented 
at  this  time  [Cardenas  and  Karplus,    1970]  by  MIMIC,    DSL/90,    and 
CSMP.     DSL/90  and  CSMP,    which  is  based  on  DSL/90,    have  transla- 
tors which  convert  the  prescribed  input  language  into  FORTRAN  sub- 
routines.    These  subroutines  are  then  compiled  by  the  FORTRAN 
compiler  and  executed  with  the  solution  routine.     This  is  a  multi-step 
process  which  is  both  time  consuming  and  not  compatible  with  the 
interactive  approach.     CSMP  and  DSL/90  will  accept  the  problem 
description  either  in  the  block  diagram  form  or  the  differential  equa- 
tion form.     There  is  a  version  of  CSMP  (similar  to  the  above  CSMP 
in  name  only)  which  uses  the  block  diagram  approach  only  and  is 
suitable  for  use  on  the  small  computers  of  the  IBM  line.     MIMIC 


incorporates  its  own  compiler  and  has  a  very  flexible  algebraic 
capability. 

The  program  described  in  this  thesis  uses  the  differential  equation 
approach  and  has  a  built-in  compiler  flexibility  similar  to  MIMIC.     In 
addition  it  has  been  designed  from  its  conception  to  be  fully  interactive 
using  graphics  as  the  interface  between  the  user  and  the  program.     All 
problem  description  is  handled  as  data.     Thus  only  the  compiled  version 
of  the  analysis  program  is  needed  by  the  user.     This  is  accomplished 
by  the  use  of  a  special  differential  equation  description  language.     An 
automatic  compiler  (transparent  to  the  user)  built  into  the  program 
provides  the  interface  between  the  description  language  and  the  fourth 
order  Runge-Kutta  solver.     This  use  of  a  special  differential  equation 
description  language  and  the  interactive  graphics  approach  give  the 
program  considerable  flexibility  in  the  type  of  differential  equations 
it  can  handle. 

The  class  of  piece-wise  continuous  non-linear  ordinary  differential 
equations  handled  by  this  program  is  restricted  by  the  following  con- 
siderations: 

1.  The  equation  must  have  only  one  independent  variable. 

2.  The  equation  must  be  normalized  so  that  the  coefficient  and 
exponent  of  the  highest  order  derivative  are  both  one. 

3.  At  each  given  moment  during  the  solution  one  and  only  one 
differential  equation  can  be  applicable.     The  applicable  differential 
equation  must  be  well  defined  for  the  conditions  at  that  point.     Note 
that  this  allows  piece -wise  continuous  differential  equations  with 
discontinuities  but  excludes  simultaneous  differential  equations. 


These  are  the  only  restrictions  inherent  in  the  method  of  solution. 
The  implementation  described  in  this  thesis  has  been  arbitrarily 
limited  to  differential  equations  of  eighth  order  or  less.      For  a  given 
problem  there  may  be  added  restrictions  in  defining  the  differential 
equations  due  to  the  limitations  of  the  program  language  used  to  define 
the  systems  and /or  the  ingenuity  of  the  user.     However,    if  the  user 
can  define  his  system  using  the  defined  programming  language  and  the 
above  three  requirements  are  not  violated  then  the  program  will  solve 
the  system. 

Although  the  language  used  in  the  system  description  is  quite 
natural  to  most  users,    it  is  important  to  completely  specify  it.     Thus 
Chapter  III  deals  with  a  description  of  the  elements  of  the  language 
and  the  principles  of  its  use. 

Other  features  of  the  program  are: 

1.  The  program  has  the  generality  required  to  be  useful  for  a 
wide  variety  of  applications.     In  addition,    the  program  does  not  re- 
quire sophisticated  computer  techniques  on  the  part  of  the  user. 

2.  The  outputs  of  the  program  are  a  smooth  x  vs  x'  plot  (phase 
plane  in  the  second  order  case)  of  the  solution  and  printer  plots  of 
the  time  response  of  each  state  variable. 

3.  The  "window"  represented  by  the  x  vs  x'  plane  can  be  adjusted 
in  size,  scale,  and  position.  These  three  parameters  can  be  adjusted 
in  the  x  and  x'  directions  independently;  thus  the  user  has  full  flexibil- 
ity to  adjust  the  window  as  desired. 

4.  The  solution  is  accomplished  using  the  state  variable  approach 
with  fourth  order  Runge-Kutta  integration.     The  parameters  of  the 
solution  are  completely  under  the  control  of  the  user. 
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5.  For  second  order  time  independent  systems  there  is  super- 
imposed on  the  phase  plane  a  grid  of  phase  trajectory  segments. 
These  slope  markers  represent  the  slope  of  the  trajectory  if  it  were 
to  pass  through  that  area  of  the  phase  plane.     Thus  without  obtaining 
a  solution  one  can  visually  obtain  a  feel  for  what  the  system  will  do 
for  various  initial  conditions.     This  in  turn  aids  the  user  in  deciding 
what  the  essential  characteristics  of  the  system  are. 

6.  For  the  Interactive  Graphics  Version,    the  parameters  of  each 
feature  mentioned  above  are  available  for  modification  as  desired  on 
an  interactive  basis.      For  example  one  may  select  specific  values  of 
various  parameters  and  observe  their  effect  on  the  resultant  solution. 
The  program  is  interactive  and  prompts  the  user  in  the  choice  of 
responses  indicating  his  desires. 

There  are  two  versions  of  this  program.      One  version  is  called 
the  Batch  Version  and  runs  on  the  IBM  360/67.      The  other  is  an 
Interactive  Graphics  Version  and  is  implemented  on  the  Electrical 
Engineering  Department's  XDS  9300  with  the  associated  ADAGE 
Graphics  Terminals.     The  interactive  graphics  version  was  primarily 
practical  due  to  the  power  of  the  "intelligent"  Adage  Terminals. 

In  Chapter  II  a  general  description  of  the  program  and  its  various 
parts  are  discussed.     A  major  portion  of  the  coding  is  the  compiler 
and  interpreter  which  are  one  of  the  primary  features  of  this  program. 

Chapter  IV  concerns  itself  with  a  description  of  the  program  from 
the  users  viewpoint.      In  addition  to  the  general  comments,    there  are 
included  sections  which  deal  separately  with  the  Batch  Version  and 
the  Interactive  Graphics  Version. 
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The  use  of  the  Interactive  Graphics  Version  requires  a  knowledge 
of  how  to  use  the  graphics  system.     This  is  explained  by  a  suggested 
first  session  exercise.     In  this  exercise  the  user  is  led  step-by-step 
through  all  the  techniques  required  to  use  the  program. 

Chapter  V  ties  all  the  above  together  by  presenting  a  large  number 
of  actual  problems  that  have  been  studied  with  this  program.     The 
input  deck  is  explicitly  presented  along  with  the  resultant  output. 

A  special  benefit  of  the  analysis  feature  of  this  program  is  its 
usefulness  as  a  teaching  or  laboratory  tool.     This  is  especially  true 
for  the  Interactive  Graphics  Version  used  to  analyse  a  second  order 
time  independent  system.      For  example  the  user  could  type  in  the 
description  for  a  simple  second  order  pendulum  system.     The  slope 
lines  give  a  very  dramatic  overall  view  of  the  effect  of  various  types 
and  amounts  of  damping.     Since  the  system  is  interactive,    the  user 
can  immediately  pursue  and  verify  theoretical  observations. 
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•m  U.       GENERAL  DESCRIPTION  OF  THE  PROGRAM 

A.      INTRODUCTION 

The  program  consists  of  five  primary  parts.     The  first  four  are 
the  compiler,    interpreter,    slope  line  generator  (used  in  the  second 
order  time  independent  case  only),    and  the  solution  generator.      Each 
of  these  will  be  explained  in  more  detail  later  in  this  chapter.      For  the 
present,    a  brief  functional  description  will  be  given. 

The  "Compiler"  takes  the  differential  equation  or  system  of  dif- 
ferential equations,    as  the  case  may  be,    and  through  a  process  to  be 
described  later  comes  up  with  a  polish  string  [Forsythe,    1969,    Bur- 
roughs,   1964]  of  integer  codes  stored  as  a  single  vector.     This  polish 
string  is  the  compiled  version  of  the  input  system  of  equations. 

The  "Interpreter,  "  when  directed,    uses  the  polish  string  and  the 
current  value  of  all  its  required  arguments  and  produces  a  value  for 
the  derivative  of  each  state  of  the  system  with  respect  to  time.      The 
required  arguments  are  indicated  by  the  system  equations  and  may 
include  the  current  values  of  the  states,    time,    and  the  values  assigned 
the  constant  coefficients  A  through  H. 

The  resultant  state  variable  derivatives  as  defined  by  the  inter- 
preter are  needed  by  both  the  slope  line  generator  and  the  solution 
generator.     The  slope  line  generator  calculates  the  grid  of  slope  lines 
as  shown  in  Chapter  V.      The  solution  generator,    as  the  name  implies, 
generates  the  time  solution. 

The  fifth  primary  part  of  the  program  is  the  coordinator.      It  is 
the  coordinator  that  ties  the  above  together  and  makes  it  work.     It 
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provides  the  required  scheduling  based  on  the  results  of  the  user's 
modifications.      These  modifications  are  accomplished  on  an  interactive 
basis  when  using  the  XDS/AGT  Interactive  Graphics  Version  and  via 
data  cards  on  the  IBM  360/67  Batch  Version. 

B.      MATHEMATICAL  BASIS  OF  THE  PROGRAM 
1.     The  General  Mathematical  Concept 

Assume  an  n      order  differential  equation  of  the  form 

,n 
2-2     +f(x)  =  g(x>u)  (2-1) 

dt 

where    f(x)  is  a  function  of  the  n-1  derivatives,    and  g(x,  u)  is  a  function 
of  the  inputs  and  the  n-1  derivatives.     This  equation  can  be   rewritten 
as  n  first  order  differential  equations.     Let  us  assume  the  following 
notation  and  substitutions 

x  =  Z(l) 

x'  =  Z(2)  =  ZDOT(l) 

x"  =  Z(3)  =  ZDOT(2) 

x,M  =  Z(4)  =  ZDOT(3)  (2-2) 


x(n_1)      =    Z(n)    =  ZDOT(n-l). 

Rewriting  equation  2-1  using  the  substitutions  and  notation  specified 
in  Equations  2-2  gives 

ZDOT(l)        =     Z(2) 
ZDOT(2)        =    Z(3) 

;  (2-3) 
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ZDOT(n-l)     =    Z(n) 
ZDOT(n)  =     g(x,u)  -  f(x). 

Each  of  the  Z's  will  be  called  the  state  variables  of  the  system. 
Note  that  for  the  first  n-1  state  variables  the  derivative  of  the  i       state 

variable  of  the  system  is  equal  to  the  i+1  state  variable.     The  deriva- 

■f"Vi 
tive  of  the  n      state  variable  is  equal  to  g(x,  u)  -  f(x). 

The  solution  generator  needs  only  these  derivatives  of  the  state 

variables  of  the  system.     The   compiler  and  interpreter  use  the  above 

observations  as  the  basis  for  the  method  of  implementation. 

Assume  that  equation  (2-1)  has  been  presented  to  the  system. 

dnx 
The  first  action  of  the  compiler  is  to  recognize  the  term     ,  n      which 

is  denoted  on  the  data  card  as  an  "x"  followed  by  n  apostrophes.     This 
notation  is  completely  analagous  to  the  use  of  dots  to  denote  time  de- 
rivatives.    The  notation  will  be  explained  in  detail  in  Chapter  III. 
Having  recognized  this  term  it  now  knows  the  order  of  the  differential 
equation.      Next  the  expression  g(x,  u)  -  f(x)  is  compiled.     The  compiled 
version  will  be  used  by  the  interpreter  to  evaluate  the  desired 
derivatives. 

The  interpreter  first  evaluates  the  compiled  version  of  g(x,  u)  -  f(x' 
which  becomes  the  value  of  the  highest  order  state  derivative.     The 
other  derivatives  are  defined  using  the  state  values  as  defined  in 
Eq.    (2-3). 

2.     Mathematical  Derivation  of  the  Slope  Lines 

The  slope  lines  are  used  only  for  second  order  time  indepen- 
dent systems.     Thus  the  comments  of  this  section  will  apply  to  a 
second  order  system  only.     The  general  system  can  be  written  as 
xn  +  f(x)  =  g(x)  (2-4) 
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where  f(x)  +  g(x,  u)  are  as  defined  in  equation  2-1. 

Letting  x  =  Z(l)  and  x'  =  Z(2)  gives 

dx/dt    =    ZDOT(l)  =  Z(2) 

(2-5) 
dx'/dt    =    ZDOT(2)  =  g(x,  u)  -  f(x). 

For  a  given  time  t  the  slope  of  the  solution  trajectory- 
is 

slope     =    (dx/dt)/(dx'/dt)  _    ,. 

=    ZDOT(2)/Z(2). 
Thus  knowing  a  position  and  a  corresponding  slope  the  desired  lines 
can  be  generated  using  simple  trigonometry. 

C.      COORDINATOR 

Figure  2-1  is  a  simplified  block  diagram  of  the  program  and  out- 
lines the  function  of  the  coordinator.     The  coordinator  provides  the 
overall  coordination  and  interface  for  obtaining  the  desired  chain  of 
events.     The  program  is  best  summarized  by  reviewing  this  block 
diagram. 

The  first  action  of  the  program  is  to  obtain  the  system  and  plot 
description  data.     This  includes  the  system  equations,    the  plot  para- 
meters such  as  scales  and  sizes,    and  the  solution  information  such  as 
the  initial  conditions,    initial  and  final  times.      With  this  information  a 
complete  run  through  the  program  can  be  made. 

Next  the  compiler  compiles  the  system  description  data.     The 
results  of  this  compilation,    if  successful,    are  saved  in  a  vector  (a 
1-dimensional  array  or  stack).      If  the  result  of  the  compilation  or  one 
of  the  steps  which  uses  the  compiled  version  encounters  an  illegal 
situation  then  an  appropriate  message  is  presented  and  the  program 
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goes  on  to  get  new  data.      For  the  Batch  Version  the  next  new  problem 
would  be  started.      For  the  Interactive  Version  the  user  would  make 
the  appropriate  corrections  and  proceed. 

After  a  successful  compilation  the  plot  is  started  and  the  grid  of 
slope  lines  is  computed  and  plotted.     The  time  solution  is  then  com- 
puted and  plotted. 

At  this  time,    the  program  is   ready  for  the  next  set  of  data  from 
the  user.      Depending  on  what  new  data  is  obtained  several  courses  of 
action  are  taken.       If  the  system  description  equations  are  changed 
then  a  recompilation  is  required.     The  program  flow  after  recompila- 
tion  is  the  same  as  before.     If  the  system  equations  were  not  changed, 
then  the  coordinator  checks  for  a  change  in  the  constant  coefficients 
or  the  plot  parameters.     A  change  in  these  will  cause  a  transfer  to 
the  slope  generating  segment. 

The  Batch  Version  differs  slightly  at  this  point.     Instead  of  asking 
if  the  "coefficients  or  plot  parameters  changed?  ,  "  the  question  should 
read:     "Is  the  next  solution  desired  on  a  new  graph?"     This  allows  a 
family  of  solutions  to  be  plotted  on  the  same  graph  if  so  desired. 

The  above  cycle  of  events   continue  until  the  user  is  done. 

D.      DESCRIPTION  OF  THE  COMPILER 

The  programming  details  of  the  compiler  are  not  of  central  im- 
portance to  this  thesis.     A  copy  of  the  coding  is  in  the  overall  program 
printout  included  as  Programs  A  and  B.      Frequent  use  of  comment 
statements  including  explanation  of  codes,    key  variables,    and  such 
are  included. 

Referring  to  the  simplified  block  diagram,    fig.    2-2,    a  very  brief 
description  of  the  compiler  will  be  given  here.     The  compiler  accepts 
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the  input  a  line  at  a  time.     The  recognition  of  the  elements  of  the 
language  is  accomplished  by  brute  force.      For  example,    let  us  assume 
that  a  "T"  has  been  recognized.     This  T  could  be  the  variable  T  or 
the  first  letter  of  either  "TAN"  or  "THEN".     A  check  of  the  following 
characters  quickly  reveals  which  of  the  three  it  is.      When  a  language 
element  is  recognized  a  corresponding  integer  code  is  placed  in  a 
stack.     At  the  end  of  the  recognition  phase,    the  entire  set  of  equations 
will  have  been  coded  into  a  string  of  integer  codes  stored  in  the  stack. 

A  blank  card  or  line  signals  that  the  recognition  phase  has  been 
completed.     At  this  time  the  string  is  in  the  form  known  as  "infix.  " 
The  final  task  of  the  compiler  is  to  convert  this  infix  string  to  a 
"polish"  string.     The  polish  string  has  the  advantage  that  it  is  quite 
easy  to  implement  an  interpreter  for  it.     This  interpreter  will  be 
discussed  in  the  next  section. 

A  part  of  the  compilation  work  is  to  recognize  and  use  information 
such  as  the  order  of  the  differential  equation.      The  order  of  the  dif- 
ferential equation  is  determined  from  the  equations  and  is  required  by 
the  interpreter. 

E.      DESCRIPTION  OF  THE  INTERPRETER 

The  purpose  of  the  Interpreter  is  to  take  the  compiled  code  which 
was  generated  in  the  compiler  and  values  for  the  required  operands 
which  were  specified  by  the  user  and  to  provide  the  derivatives  that 
are  needed  by  the  slope  and  solution  generators.      The  use  of  these 
derivatives  will  be  explained  in  the  next  two  sections. 

The  functional  block  diagram  is  presented  in  figure  2-3.      Upon 
entry  to  the  Interpreter  all  the  required  operands  have  been  defined. 
These  are  the  present  values  of  each  of  the  state  variables,    the  values 
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of  the  constants  (A  through  H),    and  the  problem  time  T.     The  Inter- 
preter takes  one  by  one  each  coded  entry  in  the  polish  stack  (result 
of  the  compilation  phase)  and  carries  out  the  required  operation.      The 
various  types  of  actions  are  summarized  by  the  block  diagram. 

The  code  may  be  carrying  information  such  as  the  order  of  the 
differential  equation  and  whether  the  differential  equation  to  be  used 
is  a  function  of  time.     These  codes  cause  the  proper  variables  and 
flags  to  be  set. 

The  logical  branching  is  done  in  the  following  manner.      Each 
statement  in  the  stack  is  separated  by  a  code  representing  a  ";".     If 
the  interpreter  evaluates  the  argument  of  an  IF  statement  to  be  false, 
then  it  skips  through  the  code  to  the  entry  following  the  next  ";". 
Thus  it  has  skipped  the  part  of  the  statement  indicating  action  to  be 
accomplished  if  the  logical  argument  is  true.     The  SKIP  instruction 
uses  the  same  method  but  will  skip  several  instructions  as  requested 
by  the  programmer. 

The  arithmetic  evaluation  of  the  polish  stack  is  accomplished  by 
the  use  of  software  LIFO  (last  in  first  out)  stacks.     The  basic  prin- 
ciples involved  in  evaluating  a  polish  stack  are  well  known  [Rosen  1967, 
Forsythe  1969].      For  ease  in  implementation  three  different  LIFO 
stacks  have  been  used.      One  is  used  for  the  arithmetic  operations,    one 
for  logical  operations,    and  one  to  hold  addresses.     Thus  they  are 
respectively  real,    logical,    and  integer  stacks.     Associated  with  these 
stacks  there  are  the  required  pointers  to  keep  track  of  their  top 
positions. 

When  the  interpreter  notes  that  a  differential  equation  has  been 
evaluated  it  terminates  the  interpretation  phase.     As  will  be  explained 
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in  Chapters  III  and  IV  the  program  uses  the  first  differential  equation 
whose  conditions  for  use  are  met.     There  it  will  be  shown  that  the 
program  easily  handles  a  system  whose  behavior  must  be  described 
by  different  differential  equations  in  different  portions  of  the  solution 
space.     The  last  action  of  the  Interpreter  is  to  define  the  desired 
system  derivatives  as  described  in  Section  B.     These  derivatives  plus 
the  order  of  the  system  are  passed  back  to  the  calling  slope  or  solution 
generator. 

F.      DESCRIPTION  OF  THE  SLOPE  GENERATOR 

The  Slope  Generator  presents  a  grid  of  slope  markers  over  the 
area  of  the  phase  plane.     The  use  of  these  slopes  is  demonstrated  in 
Chapter  V.     The  mathematical  basis  of  the  slopes  was  described  in 
Section  B  of  this  chapter. 

For  each  grid  point  on  the  phase  plane  the  interpreter  is  called. 
If  the  system,    for  that  grid  point,    is  not  second  order  or  time  indepen- 
dent then  it  skips  on  to  the  next  grid  point.     The  effect  of  this  skip  is 
to  suppress  the  slope  markers  at  that  point.      If  the  system  is  of  second 
order  and  not  time  dependent  then  the  slope  is  defined  and  is  drawn 
at  the  current  grid  point. 

The  resultant  derivative  information  from  the  interpreter  and  the 
position  information  are  used  to  calculate  the  desired  slope.      The 
slope  is  used  in  turn  to  compute  the  vector  that  represents  the  slope 
at  that  point.     This  derivation  is  made  somewhat  more  complicated 
by  the  scaling  options  available  to  the  user.      Thus  both  the  slopes  and 
the  vector  length  calculation  must  consider  the  horizontal  and  vertical 
scale  factors. 
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G.      DESCRIPTION  OF  THE  SOLUTION  GENERATOR 

.The  Solution  Generator  incorporates  a  very  straightforward  use 
of  the  Runge-Kutta  method.     The  heart  of  the  Runge-Kutta  solution 
is  the  function  subroutine  RKLDEQ.     This  function  subroutine  is  part 
of  the  IBM/360  library.      For  flexibility  the  actual  cards  are  included 
in  this  program. 

After  each  solution  point  is  generated  the  result  is  stored  for 
future  plots.     In  addition,    the  interactive  graphics  version  presents 
the  solution  as  it  takes  place.     The  solution  rate  for  a  typical  system 
is  about  ten  iterations  per  second  with  graphics.      Thus  the  resultant 
presentation  is  approximately  real  time  for  a  solution  time  step  of 
0.  1  seconds. 
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-      in.       GENERAL  DESCRIPTION  OF  THE  LANGUAGE 

A.  INTRODUCTION  TO  THE  LANGUAGE 

The  language  primarily  resembles  the  familiar  (to  most)  form  of 
FORTRAN.     Anyone  familiar  with  one  of  the  languages;  FORTRAN, 
Algol,    Basic,    or  PL/l  will  find  the  language  used  here  quite  easy. 
The  requirements  of  this  language  are  very  limited  and  of  a  special 
nature.     Thus,    since  we  have  a  special  purpose  language,    it  is  not 
complicated  by  the  generality  required  of  most  languages. 

B.  ELEMENTS  OF  THE  LANGUAGE 
1.     The  Character  Set 

The  following  characters  are  used  in  this  language. 

a.  The  alphabetic  characters  as  applicable. 

b.  The  numeric  characters,    0  through  9. 

c.  The  special  characters: 

• 

( 
) 

+ 

/ 

i 

blanks 

d.  Blanks  are  optional  and  may  be  used  as  desired  for 
clarity  and  neatness  except  in  the  middle  of  a  word,    variable  name, 
or  number.     The  semicolon  is  used  as  a  delimiter  allowing  more  than 
one  statement  per  line  and  may  optionally  be  added  at  the  end  of  a  line. 
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2.  Constants 

A  constant  is  a  fixed,    unvarying  quantity.     The  constants  used 
in  this  program  are  of  two  general  types.     They  can  either  be  self- 
defining  numeric  values  appearing  in  the  source  statements,    or  con- 
stants defined  externally  to  the  program. 

There  are  eight  externally  defined  constants  with  names  A, 
B,    C,    ...    H.     The  values  of  these  real  constants  are  defined  by  the 
user  externally  to  the  system  statements.     Thus  the  values  of  these 
constants  may  be  changed  as  desired,    between  executions,    without 
recompilation  of  the  systems  equations.     The  method  of  defining  these 
constants  is  explained  in  Chapter  IV. 

The  self-defined  constants  are  recognized  by  the  compiler 
during  compilation.      With  the  exception  of  integer  powers  all  self- 
defined  constants  are  treated  internally  as  real.     However,    when 
used  in  the  program,    the  use  of  a  decimal  point  is  optional.     That  is, 
the  decimal  point  may  be  either  explicit  or  implicit. 

In  the  case  of  exponentiation,    an  integer,    positive  or  negative, 
used  as  a  power  will  yield  the  expected  integer  power  results.     On  the 
other  hand,    non-integer  powers  are  computed  using  the  natural  log- 
arithm/exponentiation method. 

Examples: 
Valid  real  numbers:     9.  05,    0.  06,    78,    etc. 
Valid  use  of  integers:     9.7**4,    Vl2**-3,    V12  **  -  3,    etc. 
Valid  use  of  non-integer  exponents:     V23**0.  897,    etc. 

3.  Variables 

A  variable  is  a  named  quantity  whose  value  may  change 
during  execution  of  the  program.     Variables  are  specified  by  name. 
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The  name  must  start  with  the  letter  "V"  followed  by  one  or  two 
decimal  digits  representing  the  number  1  to  40  inclusive.     These 
variables  are  always  considered  as  real  variables. 

Example:    VI,    V23,    V39,    etc. 

In  addition  the  letter  "T"  is  used  as  a  variable.     T  is  the  real 
variable  which  represents  the  problem  time;  T  is  the  time  variable  of 
the  system  solution.      "X"  is  the  variable  of  the  differential  equation 
and  will  be  explained  later. 
4.       Expressions 

An  expression  is  a  string  consisting  of  constants,    variables, 
and  function  calls  interspersed  with  operators  and  parentheses.     The 
formulation  of  an  expression  is  governed  by  normal  mathematical 
convention  and  the   rules  given  below.     There  are  two  kinds  of  expres- 
sions; arithmetic  and  logical.     The  value  of  an  arithmetic  expression 
is  a  real  number.     The  value  of  a  logical  expression  is  always  a  truth 
value,    TRUE  or  FALSE.     Expressions  may  appear  in  If  Statements, 
Assignment  Statements  and  in  Differential  Equation  Statements. 

a.     Arithmetic  Expressions 

(1)    The  arithmetic  operators  are  as  follows: 
Arithmetic  Operator  Definition 


** 

Exponentiation 

* 

Multiply 

/ 

Division 

+ 

Addition 

_ 

Subtraction 

(2)    Rules  for  the  Construction  of  Arithmetic  Expressions: 
(a)    All  desired  computations  must  be  specified 
explicitly.     For  example,    to  multiply  the  variables  VI  and  V3  one  must 
write  V1*V2  and  not  VI V2. 
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(b)  No  two  numeric  operators  may  appear  in  sequence. 
(Note:     **  is  considered  as  a  single  operator).     Thus  A*  -  V19  is 
illegal  but  may  be  corrected  by  the  addition  of  parenthesis  obtaining 
A*(-V19). 

(c)  Parenthesis  may  be  used  as  desired  for  clarification 
and  to  force  the  proper  order  of  computations. 

(d)  A  mathematical  function,    together  with  its  argument, 
may  be  used  and  is  treated  as  a  variable. 

(3)  Order  of  computation  in  Arithmetic  Expressions: 

Operation  Hierarchy 

Evaluation  of  functions  1st 

Exponentiation  2nd 

Multiplication  and  division  3rd 

Addition  and  subtraction  4th 

(4)  Examples  of  Arithmetic  Expressions: 
(A*V3  +  5.  7)/EXP(Vl*B*T) 

-A*(l.  0  -  X**2)*X'  -  B*X 
A*X'*ABS(X')  +  SIGN(X) 

b.     Logical  Expressions 

(1)    Relational  Operators.     There  are  six  relational  operators, 

They,    and  the  logical  operators  which  follow,    must  always  start  and 

end  with  a  period.     The   relational  operators  are  binary  operators 

which  operate  on  arithmetic  expressions  and  express  an  arithmetic 

condition  which  can  be  either  true  or  false.     The  relational  operators 

are  as  follows: 

Relational  Operator  Definition 

.  GT.  Greater  than 

.  GE.  Greater  than  or  equal  to 
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.  LT.  Less  than 

.  LE.  Less  than  or  equal  to 

.  EQ.  Equal  to 

.  NE.  Not  equal  to 

Examples 

VI     .LT.     X' 

0.94    .  GE.     EXP(A«-  0.5) 

X    .LT.     -A    .OR.     X    .  GT.  A 

(2)  Logical  Operators.     The  three  logical  operators  are 
.NOT.     ,      .AND.      ,     and    .OR.     .     Their  use  is  defined  in  the 
ordinary  sense  of  logic. 

(3)  Rules  for  Construction  of  Logical  Expressions: 

(a)  Relational  operators  may  only  operate  on  expres- 
sions which,    when  evaluated,    have  a  numeric  value. 

(b)  Logical  operators  may  only  operate  on  expressions 
which,    when  evaluated,    have  a  logical  value  (i.  e.  ,    TRUE  or  FALSE) 

(c)  Relational  operators  may  not  appear  in  sequence. 

(d)  Two  logical  operators  may  appear  in  sequence  only 
if  the  second  one  is  the  logical  operator    .  NOT. 

(e)  A  relational  and  a  logical  operator  may  not  appear 
side  by  side  in  an  expression. 

(f)  As  in  Arithmetic  Expressions,    parenthesis  maybe 
added  as  desired  for  clarification  and  to  force  the  proper  order  of 
calculations. 

(g)  If  the  logical  operator    .  NOT.      operates  on  an  ex- 
pression with  more  than  one  term,    the  expression  must  be  enclosed 
with  parenthesis. 
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(4)  Order  of  Computations  in  Logical  Expressions. 

Operation  Hierarchy 

Evaluation  of  functions  1st 

Exponentiation  2nd 

Multiplication  and  division  3rd 

Addition  and  subtraction  4th 

.LT.,     .  LE.  ,     .  EQ.  ,     .NE.,     .  GT.  ,     .  GE.  5th 

.NOT.  6th 

.  AND.  7th 

.OR.  8th 

If  two  operators  have  the  same  hierarchy,    then  they  are  used  in  the 

order  they  appear  left  to  right. 

(5)  Examples  of  Logical  Expressions. 
(Vl**2.3     .  GE.     A)    .AND.     (X1     .  LT.      0) 

T    .GT.     1.5 

T    .GT.     1.5    .OR.     (X    .NE.     0.5)    .AND.     .NOT.     (X    .  LT.     X') 
(T     .GT.      1.5     .OR.      (X    .  NE.      0.5))    .AND.      .NOT.     (X    .  LT.     X') 
These  last  two  are  not  equivalent  since  the  added  set  of  parenthesis 
has  altered  the  computation  order  in  the  last  case. 
5.      Mathematical  Functions 

The  language  provides  for  the  use  of  various  mathematical 
functions.      The  present  library  is  listed  below.     One  calls  these 
functions  in  an  implicit  manner  by  simply  using  the  proper  function 
name,    followed  by  its  argument  in  parenthesis,    as  an  ordinary 
variable.     The  argument  itself  may  be  any  valid  mathematical  ex- 
pression which  in  turn  may  contain  function  calls. 

a.     Examples  of  Use. 

EXP(Vl) 

EXP(V1*V3/SQRT(SIN(X/X'))) 
(SIN(X)  +  COS(X**2))/3.  75 
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b.     The  various  available  functions  and  their  preas signed 

names  are  expressed  below. 

Name  Description  of  Function 

SIN  Sine  of  an  angle  in  radians 

COS  Cosine  of  an  angle  in  radians 

TAN  Tangent  of  an  angle  in  radians 

ABS  The  absolute  value  of  the  argument 

EXP  Exponential  of  the  argument 

LN  The  natural  logarithm  of  the  argument 

LOG  The  base  10  logarithm  of  the  argument 

INT  Truncates  the  real  number  to  an  integer  then  reconverts 

it  to  a  real. 

SIGN  Magnitude  of  1.  0  with  the  sign  of  the  argument. 

6.     Arithmetic  Assignment  Statements 

The  goal  in  developing  the  above  language  elements  was  to 

make  use  of  them  as  the  building  blocks  for  Arithmetic  Assignment 

Statements,    Control  Statements,    and  Differential  Equation  Statements, 

This  section  will  discuss  the  Arithmetic  Assignment  Statement. 

The  end  result  of  using  an  Arithmetic  Assignment  Statement 

is  to  compute  a  new  value  for  a  variable.      The  general  form  of  an 

Arithmetic  Assignment  Statement  is  a  =  b 

where:      a  is  a  valid  variable  name  (written  without  a  sign). 

b  is  any  mathematical  expression  following  the  above 

guidelines. 

=  is  an  order  to  compute  the  value  of  the  expression  b 

on  the  right,    and  then  give  that  value  to  the  variable 

named  on  the  left. 

A  lengthy  discussion  can  be  made  on  the  difference  between  the 

Arithmetic  Assignment  Statement  and  a  mathematical  statement  of 
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equality.     It  is  perhaps  unfortunate  that  an  "="  is  traditionally  used  in 
both. cases.     However  if  one  keeps  the  above  definition  in  mind,    there 
should  be  no  confusion.      For  someone  completely  new  to  computer 
programming  methods,    it  is  recommended  that  one  of  the  many  good 
fundamental  FORTRAN  manuals  be  consulted.     The  formation  of  and 
use  of  the  Arithmetic  Assignment  Statement  is  identical  with  that  of 
FORTRAN. 

Examples: 

1  2 

Let  VI  =~    + 


2x  .    2 

4x 


then  the  coded  version  is: 

VI  =  -1.  /(2.  *X)  +  A**2/(4#X**2) 

,  2/5 
Let  V5    =    (A  +  Bx    ) 

then  the  coded  version  is: 

V5  =  (A  +  B*X**3)**(2.  0/5.  0) 

?  2     ?- 

Let  V9  =  (SIN  x  +  COS   x)2 

then  the  coded  version  would  be: 

V9  =  (SIN(X)**2  +  COS(X)**2)**0.  5. 


Let  V23  = - —  +  log       tan   ^ 

cosx  °     I  L    \ 

then  the  coded  version  would  be: 

V23  =  1.  0/COS(X)  +  LOG(ABS(TAN(X/2.  0))). 

Let  V5  =  V5   .  ,    +   lax2  I 

new  old         J  1 

then  a  coded  version  would  be: 

V5    =    V5  +  ABS(A*X*X). 
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7.      Control  Statements 

a.     The  Skip  Statement 

The  control  statements  are  of  two  general  types.     One 
is  a  simple  unconditional  transfer  statement  called  a  Skip  Statement. 
Its  form  is: 

SKIP    N 
where:     N  is  an  integer  that  specifies  how  many  statements   should  be 
skipped. 

b.  The  Arithmetic  IF  Statement 

The  other  more  powerful  control  statement  is  the 
Arithmetic  IF  Statement.     This  is  essentially  the  statement  that  one 
is  familiar  with  in  FORTRAN.     Its  form  is: 

IF(a)    THEN  b; 
where:    a  is  any  valid  logical  expression  as  previously  defined, 
b  may  be  an  Arithmetic  Assignment  Statement,    a  Skip 
Statement,    or  a  Differential  Equation  Statement. 
It  commands  the  computer  to  evaluate  the  logical  expression  a.     If  a 
is  true  then  do  b,    otherwise  skip  b  and  go  on  to  the  next  statement. 
The  use  of  the  keyword  THEN  is  optional  as  is  the  ";". 

c.  Examples  of  Valid  Control  Statements 

SKIP  3 

IF(X    .GT.  X'   )    THEN  SKIP  4 

IF(X    .  LE.  ABS(X'))    THENV1=A*X«=2 

IF(X    .  GE.  0.0)    THEN    X"  +  0.  2*X'     +    X    =0.0 

8.     Differential  Equation  Statements 
a.      Notation 

Differential  Equation  Statements  are  intended  to  be  as 
natural  as  possible.     The  apostrophe  notation,    which  is  widely  used, 
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has  been  adopted.     Thus  this  program  used  x"    as  the  symbol  for  the 
second  derivative  of  x. 

Other  derivatives  of  X  are  handled  in  an  analogous  manner. 
With  the  exception  of  this  notation  for  the  derivatives,    the  remainder 
of  the  differential  equation  is  written  using  the  construction  methods 
used  for  mathematical  expressions.      Note  that  in  the  case  of  a  Dif- 
ferential Equation  Statement,    the  statement  actually  expresses  an 
equation.     Thus  for  this  statement  only,    the  equal  sign  denotes 
equality  rather  than  an  assignment. 
Example:    X,M  +  A*X'*ABS(X')  +  X**2  =  V1*T 

b.  Writing  a  Differential  Equation  Statement 

The  form  of  a  Differential  Equation  Statement  is: 

a  —  b  =  c 
where:      a  denotes  the  highest  derivative  of  X  written  using 
the  notation  of  the  previous   section, 
b  and  c  denote  any  valid  mathematical  expressions 
using  any  or  all  of  the  constants  and  variables  and 
using  any  or  all  of  the  lower  derivatives  of  X  as 
variables. 
It  can  be  seen  that  the  rules  for  writing  b  and  c  are  extremely  liberal. 
Thus  the  type  of  differential  equation  that  can  be  handled  is  very 
general.      The  differential  equations  may  be  extremely  non-linear, 
and  have  coefficients  that  vary  with  time  and/or  position. 

c.  Examples 

The  above  generality  can  best  be  demonstrated  by  a  few 
examples.     Valid  constructions  are: 
X'"    +    A*X"     +    0.75*X'    +    0.0 
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X"    +    (X*X')**X    +    SIN(X)    =  10.0 

x..n    +    x.n    +    X'7X    =    LOG(X"'**0.  5) 

C.      USING  THE  LANGUAGE 

1.  General  Comments 

The  goal  of  all  the  above  elements  of  the  language  has  been  to 
develop  tools  for  use  in  defining  the  differential  equations  of  a  system. 
The  above  language  features  make  it  possible  to  handle  a  wide  variety 
of  system  types.     This  section  demonstrates,    by  use  of  several  ex- 
amples,   how  to  use  the  elements  of  the  language  to  describe  various 
systems.     The  series  of  examples  will  start  with  simple  applications 
and  then  progress  into  more  involved  descriptions. 

2.  A  Simple  Differential  Equation 

The  simplest,    but  still  very  useful  and  powerful,    form  of 
system  description  consists  of  one  Differential  Equation  Statement. 
Let  us  take  for  example  Van  der  Pol's 
Equation:  X"  -  A*(l.  0  -  X**2)*X'    +  B*X  =  0.  0 

This  one  data  card,    or  line  on  the  graphics  terminal,    completely 
describes  the  differential  equation  with  the  exception  of  A  and  B. 
The  variables  A  and  B  are  specified  on  an  additional  data  card.      When 
it  is  desired  to  alter  the  values  of  A  and  B  they  are  changed  independ- 
ently of  the  system  equations.      Considerable  execution  time  is  saved 
by  eliminating  re  compilation  of  the  system  equations  when  it  is  only 
desired  to  alter  the  value  of  a  coefficient.     If,    however,    one  did  not 
anticipate  varying  coefficients  A  and  B,    then  the  differential  equation 
could  be  coded  as  follows,    for  example: 
X"     -    0.  95*(1.  0  -  X**2)*X'    +    1.7*X    =0.0 
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Assume  now  that  it  was  desired  that  the  first  coefficient  be  0.  90  vice 
0.95.      There  is  still  no  problem,    the  required  recompilation  is  trans- 
parent to  the  user.     The  on  line  user  will  only  notice  a  somewhat 
longer  delay  time  before  the  solution  starts. 

3.     More  than  One  Differential  Equation  Statement 

In  most  control  system  applications  more  than  one  differential 
equation  is  needed  to  describe  a  given  system.     However,    only  one 
differential  equation  is  applicable  at  a  given  moment  in  the  solution. 
For  example,    due  to  a  particular  type  of  non-linearity,    the  selection 
of  the  proper  differential  equation  may  be  dependent  on  the  sign  of  the 
velocity  or  position.      In  another  situation,    the  selection  of  the  dif- 
ferential equation  may  be  a  function  of  time.      Numerous  examples  of 
various  control  system  descriptions  appear  in  Chapter  V. 

For  purposes  of  illustration,    we  shall  assume  the  following 
system: 

If  X  is  less  than  -0.  2 
then  the  differential  equation  is 

X"     +    0.  2*X'     +  X     =0.2 

If  X  is  greater  than  0.  2 
then  the  differential  equation  is 

X"    +    0.  2*X'    +    X    =  -0. 2 

If  X  is  neither  of  the  above 
then  the  differential  equation  is 
X"    +    0. 2*X'    +    X    =  0 

The  above  system  could  be  coded  as 

IF(X    .  LT.     -     0.2)    THEN  X"    +    0.  2*X'     +    X    =0.2 
IF(X    .GT.     0.2)    THENX"    +    0.2*X'    +    X    =-0.2 
X"    +    0.  2*X'    +    X    =  0.  0 
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During  the  execution  of  the  solution,    when  the  differential 
equation  is  referred  to,    the  Interpreter  will  test  the  argument  of  each 
IF  statement  starting  at  the  top.      It  will  use  the  differential  equation 
corresponding  to  the  first  condition  it  finds  to  be  true.      Note  that  the 
third  statement  is  not  preceded  by  a  condition.     The  effect  of  the  third 
statement  is  that  the  condition  for  its  use  is  always  true.     Thus  if  the 
Interpreter  does  not  use  one  of  the  first  two  differential  equations  it 
will  use  the  third. 

To  further  clarify  the  selection  of  the  right  differential 

equation,    consider  the  following: 

X"    +    0.2*X'     +    X    =     0.  0 

IF(X    .LT.     -0.2)    THENX"     +    0.  2*X'    +    X    =0.2 

IF(X    .  GT.     0.2)    THENX"    +    0.2*X'    +    X    =    -0.  2 

This  is  the  same  set  of  statements  as  before  except  that  the 
bottom  statement  has  been  moved  to  the  top.     However,    the   result  is 
drastically  different.     The  condition  for  use  of  the  present  first  dif- 
ferential equation  is  true  by  default.     Thus  it  will  be  used  regardless 
of  the  value  of  X.     Again,    the  point  to  be  stressed  is,    only  one  dif- 
ferential equation  can  be  used  at  a  given  moment;  the  one  used  will  be 
the  first  one  with  a  true  condition. 

4.      Using  the  Remaining  Features  of  the  Language 

Often  it  is  convenient  to  be  able  to  define  and  calculate  some 
variables  before  the  Differential  Equation  Statement.      Consider  the 
following  coded  system  description: 

IF((4.  0*X  +  X')    .GE.     0.0)    THEN  X"  +  0.  02*X'    +    0.3    =0.0 
IF((4.  0*X  +  X')    .LT.     0.0)    THEN  X"  +  0.  02  *X'     -     0.3    =0.0 

Using  features  of  the  language  previously  described  this  could  be 

recoded  as  follows: 
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VI  =  +  0.  3;     IF((4.  0*X  +  X*)    .  LT.     0.0)    THEN  VI  =  -VI 
X"    +    0.  02*X'    +    VI     =    0.0 

The  semicolon  allows  more  than  one  statement  per  line. 
This   recoded  version  has  the  same  effect  as  the  previous  version. 
Note  that  there  is  only  one  differential  equation.     The  VI  used  in  the 
differential  equation  is  defined  in  the  line  above  in  such  a  way  that 
the  effect  is  the  same  as  the  two  differential  equations  of  the  first 
version. 

For  this  system  an  even  simpler  coding  is  possible. 
X"  +0.02*X!     +    0.  3*SIGN(4.  0*X    +    X')    ■    0.0.     This  example 
demonstrates  the  versatility  of  the  language  by  showing  the  same 
system  coded  in  three  different  ways. 
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IV.       PRINCIPLES  AND  METHODS  OF  USE 

A.  GENERAL  COMMENTS 

The  previous   chapters  have  concentrated  on  various  specific 
aspects  of  the  program.     In  Chapter  II  the  internal  structure  of  the 
program  was  described.      Chapter  III  described  the  language  associated 
with  use  of  the  program.      This  chapter  will  endeavor  to  pull  these  to- 
gether with  specific  instructions  on  the  use  of  the  program.     In  Section 
B  those  areas  common  to  both  the  Batch  Version  and  the  Interactive 
Graphics  Version  will  be  discussed.     Sections   C  and  D  will  deal  with 
the  peculiarities  of  the  two  versions   respectively. 

B.  GENERAL  PRINCIPLES 

1.  System  Description 

One  of  the  first  tasks  for  the  user  is  to  encode  the  system 
equations  in  the  form  usable  by  the  computer.     One  of  the  primary 
features  of  the  program  is  that  this  is  done  in  a  very  natural  way. 
For  simple  systems  this  will  involve  writing  only  a  single  equation. 
For  systems  with  a  complicated  control  strategy  the  user  may  have 
to  write  what  amounts  to  a  small  program  to  completely  specify  it. 
This   general  area  was  covered  with  examples  in  Section  C  of  Chapter 
III.      Chapter  V    will  further  amplify  this  with  a  large 
number  of  complete  examples. 

2.  Use  of  the  Constants 

The  use  of  the  constants  was  described  in  Chapter  III.     How- 
ever for  the  completeness  of  this   Chapter  they  will  be  mentioned 
again.     Often  when  preparing  to  study  a  particular  system  it  may  be 
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desired  to  vary  a  parameter  over  several  values.      One  way  to  accom- 
plish this  is  to  vary  the  numerical  value  in  the  system  equations. 
However  this  amounts  to  starting  a  whole  new  problem  and  the  system 
equations  must  be  recompiled. 

An  easier  way  to  accomplish  this  parameter  variation  is  to 
use  the  constants.     These  are  the  constants  labeled  A  -  H.     When 
writing  the  system  equation  use  the  constants  in  place  of  an  actual 
numerical  value  and  define  them  separately  as  desired.      When  it  is 
desired  to  change  the  parameter,    redefine  the  associated  constants. 

Another  handy  use  of  these  constants  is  when  the  same 
number  is  to  be  used  several  places  in  the  system  description.     In 
this  use  a  constant  would  be  defined  once  and  then  used  as  often  as 
desired. 

3.       The  Phase  Plane  Window 

The  output  of  the  program  is  a  phase  plane.     The  slopes  and 
solutions  are  presented  on  this  phase  plane.      This  phase  plane  can  be 
visualized  as  a  window  in  which  the  user  has   control  over  the  place- 
ment,   size,    and  scales  to  be  used  in  its  construction.     The  controls 
are  independent  in  the  horizontal  and  vertical  directions.      Following 
are  the  control  parameters  and  their  use: 

X-SCALE         The  horizontal  scale  in  units  per  inch. 

Y-SCALE         The  vertical  scale  in  units  per  inch. 

X- CENTER     The  X  coordinate  of  the  center  of  the  window. 

Y- CENTER     The  Y  coordinate  of  the  center  of  the  window. 

X-SIZE  The  horizontal  dimension  in  inches. 

Y-SIZE  The  vertical  dimension  in  inches. 
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4.     Solution  Parameters 

The  remaining  values  to  be  specified  are  the  parameters  of 
the  solution.  These  include  the  initial  time,  time  step,  final  time, 
and  the  initial  conditions. 

In  most  cases,    the  user  will  let  the  initial  time  be  zero  and 
the  final  time  be  some  positive  value  with  an  associated  positive  time 
step.     However,    this  is  not  required.      In  fact  the  final  time  may  be 
less  than  the  initial  time.     In  this  later  case  the  time  step  would  be 
negative. 

The  program  can  handle  systems  up  to  eighth  order.     Thus 
there  are  provisions  for  specifying  up  to  eight  initial  conditions. 
This  limitation  to  eighth  order  was  arbitrarily  set  by  the  ease  of 
handling  the  initial  values  on  input.      Eight  initial  values  fit  nicely  on 
one  input  card. 

C.      SPECIFIC  DETAILS  ON  THE  BATCH  VERSION 
1.      Planning  for  the  Batch  Run 

For  the  batch  version  the  user  must  plan  the  entire  run 
beforehand.      One  may  be  uncertain  of  some  parameters  of  the 
problem  on  the  first  run.     If  this  is  the  case,    the  user  should  take 
advantage  of  the  flexibility  of  this  program  on  the  first  submission. 
That  is,    the  first  submission  should  include  combinations  of  values 
such  that  when  the  results  are  returned  one  will  know  what  values 
to  choose  for  the  next  run.      For  example,    take  the  problem  of  the 
choice  of  time  step  and  final  time.      Here  it  would  be  good  to  make 
a  guess  at  a  suitable  time  step  and  final  time.     Then   run  an  addition- 
al two  solutions  with  these  times  multiplied  by  ten  and  divided  by 
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ten.     The  added  effort  for  these  additional  solutions  would  be  trivial. 

When  the  results  come     back  it  will  then  be  quite  easy  to  decide  on 

the  best  values. 

The  same  idea  can  apply  to  other  areas  such  as  scaling.     The 

time  required  to  punch  the  added  cards  for  a  few  extra  runs  can  save 

considerable  time  lost  to  turn-around-time  if  an  added  exploratory  run 

is  required. 

2.       Data  Card  Arrangement 

Cards   1-2  Plot  title  including  your  name  (6A8) 

Format 

3  X-SCALE(Units/inch), 

X- CENTER  (Coordinate  of  X- CENTER) 

X-SIZE  (Inches) 

Y-SCALE  (As  above) 

Y- CENTER  (As  above) 

Y-SIZE  (As  above) 

These  six  items  are  placed  on  one 

card  in  the  above  order  using  a 

(6F10.  0)  Format. 

If  (X-SIZE     .  GT.     9.0)  then  the  plot  is 

rotated  CCW  on  the  paper  by  a  quarter 

turn.      Due  to  paper  size  only  one  of 

the  dimensions  may  exceed  9.  0  inches 

and  must  not  be  more  than  15.  0  inches. 

4  SIZE  of  slope  markers  (inches,  0.2 
works  nicely),  NUMBER  per  inch  (3 
recommended) 

In  (F10.4,  11)  Format 

5,  6,  7,8...        EQUATIONS  TO  BE  USED:    Following 
the  rules  of  Chapter  m  (80A1   Format). 

Then  follow  four  cards  per  solution  desired  as 
described  below: 

-First-  VALUES  of  A,  B,  C,  D,  E,  F,  G,  H  (8F10.  4) 

Format 
-Second-         Integration  info    (INITIAL  TIME,    TIME 

STEP,    and  FINAL  TIME)  (900  Solution 

steps  Max),    (3F10.4)  Format. 
-Third-  INITIAL  CONDITIONS  for  X,    X1,    X", 

...    in  (8F10.  4)  Format 
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-Fourth-  A  "WHATS  NEXT  CARD"  use  the  following 

code: 

0  -    There  are  no  requests  following. 

1  -    The  following  are  the  second  and 

third  cards  with  new  integration 
and  initial  condition  information 
only.      I  want  this  new  solution  to 
be  plotted  on  the  present  graph. 

2  -    The  following  is  another  four  cards 

with  data  for  another  solution  of  the 
present  system  on  a  new  graph. 

3  -    The  following  card  starts  a 

completely  new  problem  starting 
with  the  plot  title  cards. 

3.     Sample  Data  Deck 

Following  is  a  typical  sample  deck  of  input  cards.     The  com- 
ments enclosed  within  parenthesis   are  not  part  of  the  cards  but  are 
included  in  this  listing  for  clarification  only.     These  cards  accomplish 
the  following. 

The  system  X"     +  A*X'  +  SIN(X)  =  0  is  compiled  and  two 
solutions  are  obtained  with  A  =  0.  3.     Then  A  is   changed  to  0.  7  and 
two  more  solutions  are  run  using  the  same  time  information  and 
initial  conditions  as  with  A  =  0.  3.      Then  a  completely  new  problem 
(a  saturated  servo  system)  is  defined. 

The  input  cards  used  are: 

NELSON,    H.    G.  Pendulum  (GRAPH  TITLE  -   1st  LINE) 

X"  =  A*X'  =  SIN(X)  =  0  (GRAPH  TITLE  -  2nd  LINE) 

2.0  0.0  8.0  2.0  0.0  6.0    (WINDOW  INFO) 

0.  2  4  (SLOPE  MARKER  SIZE  AND  NO.  ) 

Xn  +  A*X'  +  SIN(X)  =  0  (SYSTEM  EQUATION) 

(END  OF  SYSTEM  EQUATIONS) 

0.  3  (VALUE  OF  A) 

0.0  0.05  25.0  (TIME  INFO) 

-7.0  3.0  (INITIAL  CONDITIONS) 

1  (WHAT'S  NEXT  SIGNAL) 
0.0              0.05           25.0  (TIME  INFO) 

-7.0  5.0  (INITIAL  CONDITIONS) 

2  (WHAT'S  NEXT  SIGNAL) 
0.  7  (VALUE  OF  A) 

0.0  0.05  25.0  (TIME  INFO) 


43 


-7.  0           3.  0 

1 

0.  0             0.  05 

25.  0 

-7.0           5.0 

3 

NELSON,    H.    G. 

0902 

(INITIAL  CONDITIONS) 
(WHAT'S  NEXT  SIGNAL) 
(TIME  INFO) 
(INITIAL  CONDITIONS) 
(WHAT'S  NEXT  SIGNAL) 
(GRAPH  TITLE   -  1st  LINE) 
SATURATED  SERVO  SYSTEM  (GRAPH  TITLE  -   2iri  LINE) 
0.3  0.0  6.0  0.2  0.0  8.0      (WINDOW  INFO) 

0.2  4  (SLOPE  MARKER  SIZE  AND  NO.  ) 

IF(X.  LT.  (1.  0-2.  0)*0.2  THEN  X"  +  0.  2*X'  +  X  =  0.  2 
IF(X.  GT.  0.  2)  THEN  X"  +  0.2*X'  +  X  =  -0.2    (SYSTEM  EQUATIONS) 
X"  +  0.  2*X'  +  X  =  0.  0 

(END  OF  SYSTEM  EQUATIONS) 
0.0  0.05  30.0  (TIME  INFO) 

0.7  (INITIAL  CONDTT IONS) 

(WHAT'S  NEXT  =  BLANK  =  0,    THUS 
END  OF  DATA) 


D.      SPECIFIC  DETAILS  ON  THE  INTERACTIVE  GRAPHICS  VERSION 

1.  Planning  for  the  Interactive  Session 

Since  the  Interactive  Version  is  very  easy  to  use  and,    as  the 
name  implies,    is  interactive  the  required  preparation  is  minimal. 
This  statement  assumes  the  first-time  user  has   read  this  manual.      If 
a  mistake  is  made  or  a  bad  guess  is  made  for  a  parameter  value,    it 
can  normally  be  corrected  before  the  run  starts.      If  one  notices  the 
error  after  the  run  has  started,    there  is  still  no  problem.     The  pro- 
gram has  a  large  amount  of  protection  built  into  it.     Hopefully,    no 
mistakes  or  errors  by  the  user  can  cause  the  interaction  to  terminate. 
Regardless  of  the  immediate  results  of  a  mistake,    control  will  return 
to  the  user  for  desired  corrections  and/or  modifications. 

2.  Check-off  List  for  Loading  the  Program 

Following  is  a  step  by  step  list  to  be  used  to  load  and  run  the 
program. 

a.     Turn  on  the  XDS  9300  Computer  (if  not  already  on)  using 
the  instructions  attached  to  the  console.     Also  turn  on  Sense  Switch  2. 
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b.  Mount  the  binary  tape  of  the  program  on  one  of  the  tape 
drives  and  advance  it  to  the  load  point.     Set  the  "unit  select"  to  0,    the 
"density  select"  to  556  and  the  mode  switch  to  automatic.     The  "unit 
power",    "unit  ready",    "file  protect",    and  "load  point"  lights  should  be 
on. 

c.  Place  the   small  card  deck  as  specified  in  Appendix  B  in 
the  card  reader  and  press  the  "power  on"  and  "start"  buttons  in  that 
order. 

d.  Ready  the  printer.     On  the  XDS  9300,    verify  that  "Sense 

Switch  2"  is  on.     Then  press  the  following  buttons: 

HALT 

RESET 

CLEAR  and  CLEAR  FLAGS  simultaneously 

RUN 

CARDS. 

e.  The  card  reader  should  begin  reading  cards  and  then  the 
Tape  Drive  will  begin  operating.  It  will  take  about  three  minutes  to 
load  and  link  the  program. 

f.  Following  the  instructions  at  the  AGT's,  turn  on  the  AGT 
unit  to  be  used  and  load  "Gated". 

g.  Verify  that  "Gated"  is  ready  by  pressing  the  "Text  Edit" 
button  on  the  moveable  button  panel.  The  following  should  appear  at 
the  bottom  of  the  screen: 

"TEXT  BLOCK  SELECT  MODE  Block  1" 

h.     When  the  XDS  console  typewriter  types   "ENTER 
IDEV  =  1  or  2"  do  the  following: 

1.  Ensure  the  AGT  is  ready 

2.  Note  the  number  of  the  AGT 

3.  Type  "IDEV  =  1*",    return 

or  "IDEV  =  2*",    return  as  applicable. 
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i.     This  completes  the  preparation  phase.     Although  this  list 
appears  long  or  complicated  it  is  quite  easy  once  it  has  been  done  a 
few  times.     Quite  often  the  user  will  find  both  the  XDS  and  the  AGT 
already  turned  on  and  ready  to  go.     It  is  strongly  recommended  that 
the  initial  use  of  this  program  be  made  during  working  hours  when 
advice  and/or  help  is  available. 

3.       Interacting  with  the  Program 

The  interaction  starts  with  a  question  to  the  user.      "Do  you 
desire  to  start  a  completely  new  problem?,    answer  'YES'  or  'NO'  ". 
Answer  "YES"  or  "NO"  on  the  keyboard  and  then  hit  the  return  key. 
If  you  answer  "NO",    a  problem  already  in  the  computer  will  be  dis- 
played.    This  is  useful  for  becoming  familiar  with  the  program.     It 
is  suggested  for  the  first  time  user  that  the  answer  should  be  "NO" 
to  that  question.     Then  the  new  user  should  try  various  manipulations 
on  the  given  problem  until  a  working  knowledge  of  the  program  is 
obtained. 

After  each  run  the  control  returns  again  to  this  same  question. 
After  the  first  pass  the  answer  to  the  question  gives  the  user  the 
option  of  making  modifications  to  the  present  system  already  on  the 
screen  or  starting  a  complete  new  problem.      If  there  are  only  a  few 
changes  to  be  made  it  is  best  to  selectively  edit  the  present  problem. 
If  it  is  desired  to  start  a  new  problem  or  make  many  changes,    it  will 
probably  be  faster  to  use  the  "INPUT"  mode.      When  using  the  "INPUT" 
mode,    the  system  prompts  the  user  for  the  desired  information. 

The  easiest  way  to  learn  the  system  is  to  use  it.     Thus  the 
following  section  is  a  step-by-step  procedure  for  using  the  program 
on  several  specific  problems.      It  is   suggested  that  the  user  closely 
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follow  this  guide  using  the  given  problems.     Once  the  user  is  familiar 
with  the  procedures  for  using  the  program  and  the  graphics  terminal 
then  he  can  proceed  on  his  own. 

4.     Suggested  First  Session  Example 

a.  Load  the  program  as  described  in  Section  2.     The  question 
"DO  YOU  DESIRE  TO  START  A  COMPLETELY  NEW  PROBLEM?  , 
ANSWER  YES  OR  NO"  will  appear  on  the  screen. 

b.  Type  "NO",    return.     At  this  time  the  internally  stored 
program  will  appear  on  the  screen.     In  addition  will  be  the  message 
"MAKE  DESIRED  CHANGES  AND  CHOOSE  NEW  INITIAL  CONDITIONS 
(TEXT  OR  LIGHT  PEN)  TO  INITIATE  NEW  RUN".     The  signal  to  the 
program  that  the  user  is  ready  to  proceed  is  to  edit  the  first  line  of 
the  initial  conditions  or  to  select  initial  conditions   graphically  with 
light  pen. 

c.  If  satisfied  with  the  initial  conditions  already  on  the 
screen,    make  the  following  null  edit.      Press  the  following  buttons  in 
succession: 

TEXT  EDIT 

NEXT  BLOCK     -    Check  to  see  that  the  second  line  characters  are 

slanted.     If  not  use  the  NEXT  BLOCK  and/or 
PREVIOUS  BLOCK  until  line  two  is   slanted.      The 
slanting  indicates  that  line  two  is   ready  to  be 
edited. 

GO  EDIT  -    Corrections   can  now  be  made  as  desired  to  line 

two.      For  this  time  a  null  correction  will  be 
made,    i.  e.  ,    no  correction. 

TERMINATE  EDIT 

The  program  now  has  all  the  information  it  needs  to  make  a  run  and 
the  above  edit  of  line  two  has  indicated  that  the  program  is  ready  to 
be  run.  When  the  run  has  terminated  the  question  of  step  a.  will  be 
repeated.  A  complete  loop  in  the  program  has  now  been  completed, 
d.     Again  answer  by  typing  "NO",    return. 
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e.  Before  this  next  run  the  goal  will  be  to  change  the  initial 
conditions  to  -3.  8  and  3.  0  respectively.     To  do  this  "select"  and  "go 
edit"  on  line  two  as  in  step  c  above.     This  time  a  correction  will  be 
made.     While  holding  the  typewriter  key  marked  "CTRL",    space  over 
with  the  "F"  (forward)  key  until  the  underscore  is  immediately  to  the 
right  of  the  first  initial  condition  to  be  corrected.      Now  release  the 
"CTRL"  key.      Press  the  "RUB  OUT"  key  the  desired  number  of  times 
to  eliminate  the  portion  to  be  changed.     Then  type  in  the  new  first 
initial  condition.     Be  sure  to  replace  as  many  characters  as  were 
removed.     This  keeps  the  whole  line  in  the  proper  alignment.      Next, 
space  over  and  correct  the   second  initial  condition  using  the  same  keys 
as  used  to  correct  the  first  initial  condition.     When  you  are  satisfied 
with  the  line  press  the  terminate  edit  button  or  the  "RETURN"  key  on 
the  typewriter.     Either  one  will  terminate  the  edit  and  initiate  the  next 
run.     The  end  of  this   run  will  bring  you  back  to  the  question  of  step  a. 

f.  Again  answer  by  typing  "NO",    return. 

g.  The  goal  of  this  run  is  to  select  a  set  of  initial  conditions 
using  the  light  pen.      Press  in  sequence  the  "graphics  edit",    "go  edit", 
and  "cursor  on"  buttons.      Insure  that  the  light  pen  is  turned  on  by 
checking  for  the  pilot  light  on  the  box  at  the  end  of  the  light  pen  cable. 
Take  the  light  pen,    point  it  at  the  cursor,    press  the  little  white  button 
on  the  light  pen  and  lead  the  cursor  around  the  screen.      When  you  are 
familiar  with  positioning  the  cursor,    select  with  the  cursor  a  desired 
initial  condition,    then  press  the  "move"  and  then  "draw"  buttons. 
With  the  light  pen  move  the  cursor  to  another  initial  condition  on  the 
screen.     There  should  be  a  dot  remaining  at  the  location  of  the  first 
initial  condition  pair  selected.      Pressing  the  "move",    "draw"  buttons 
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will  select  the  present  cursor  position  as  a  second  initial  condition 
pair.     This  can  be  repeated  up  to  fifteen  times  if  desired.      When  you 
are  through  selecting  initial  conditions,    terminate  the  edit.     This  will, 
as  when  selecting  the  initial  conditions  numerically,    initiate  the  run. 
A  solution  will  now  be  obtained  for  each  of  the  selected  initial  con- 
ditions and  then  control  will  return  again  to  the  question  of  step  a. 

h.     Again  answer  "NO". 

i.     This  time  change  A  to  3.  0  .     To  do  this  press  "text  edit" 
and  then  press  next  block  until  the  characters  in  the  line  containing 
the  A  tilt.      You  are  now  ready  to  edit  this  line  using  the  method  of 
step  e  above. 

j.      When  through  editing  A  initiate  the  next  run  using  the 
technique  of  steps  c,    e,    or  g  above. 

k.     Once  more  answer  "NO". 

1.     This  time  it  is  desired  to  edit  the  system  equation  to 
obtain  X"  +  X*X'  +  X  =  0.      Press  the  "text  edit"  and  then  the  "next 
block"  button  about  eight  times  until  the  equation  characters  slant 
indicating  they  have  been  selected.      Make  the  modification  as  before. 
You  do  not  need  to  replace  as  many  characters  as  you  removed  when 
editing  the  equations.      When  ready  initiate  the  next  run  using  one  of 
the  previously  discussed  methods. 

m.      When  the  solution  has  been  completed  and  the  question 
of  step  a  returns  answer  "YES".     In  response  to  each  question  make 
responses  as  suggested  below.     After  each  response  hit  the  return 
key. 
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SYSTEM  EQUATIONS: 

VI  =  4.  0*X  +  X' 

IF(V1.  GE.  0.  0)       THEN  X"  +  0.  02*X'  +  0.  3  =  0.  0 

X"  +  0.  02*X"  -  0.  3  =  0.  0 

blank  line 

X-SCALE  =0.2 

X-CENTER  =0.0 

X-SIZE  =6.0 

Y-SCALE  =0.2 

Y-CENTER  =0.0 

Y-SIZE  =8.0 

A  -  H  hit  return,    will  default  to  zero 

INITIAL  TIME  =     0.  0 

TIME  STEP  =    0.  04 

FINAL  TIME  =     10.  0 

FIRST  I.    C.  =0.6 

remainder  of  I.  C.  's,    hit  return, 

n.      Now  the  program  will  give  the  question  "MAKE  DESIRED 
CHANGES  (TEXT  OR  LIGHT  PEN)  TO  INITIATE  NEW  RUN".     At 
this  point  you  can  proceed  as  in  steps  c.    through  1.    above  as  desired. 
Thus  if  any  corrections  need  to  be  made  as  a  result  of  incorrect 
responses  to  the  above  questions  they  can  be  made  at  this  time.     When 
ready  to  initiate  the  run  use  the  method  of  step  c.     It  is  best  to  select 
the  initial  conditions  numerically  for  the  first  run  after  changing  the 
window  parameters.     Initial  conditions   selected  with  the  light  pen  take 
on  the  corresponding  values  of  that  point  relative  to  the  current 
picture  on  the  screen. 

o.      You  have  at  this  time  exercised  all  the  control  functions 
needed  for  full  interactive  control  of  the  program. 
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V.       APPLICATIONS  TO  VARIOUS  PROBLEMS 

A.      GENERAL  COMMENTS 

The  purpose  of  this  chapter  is  to  demonstrate  the  use  of  the 
program  by  using  it  to  study  several  classes  of  control  system  prob- 
lems.    It  is  not  the  intention  of  this   chapter  to  derive  the  applicable 
equations,    but  rather  to  show  the  use  of  the  program  in  their  solution. 
If  the  purpose  of  the  system  is  to  reduce  the  error  to  zero  or  a 
minimum,    then  it  is  convenient  to  use  error  as  the  dependent  variable 
instead  of  the  regular  output  variable. 

In  the  following  examples  the  problem  will  be  defined  first  then 
the  data  cards  will  be  shown,    followed  by  the  resultant  phase  plane. 
Although  these  examples  will  have  been  produced  using  the  hard  plots 
of  the  Batch  Version,    the  input  data  is  identical  with  that  of  the 
Interactive  Graphics  Version. 
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B.      VARIOUS  SINGLE  EQUATION  EXAMPLES 
Example  1 
Van  der  Pol's  Equation 

X"   -  A*(l.  0  -  X**2)*X'  +  B*X  =  0.  0 

WHERE: 

A  =  1.  0 
B  =  1.  0 

PLOT  PARAMETERS: 

XSCALE  =  1.0 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes     =     0.2 

Number  of  slopes  per  inch    =    4 


YSCALE 

~ 

1.  0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

8.  0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

X"  -  A*(l.  0  -  X**2)*X'  +  B*X  =  0.  0 
1.0  0.0  6.0  1.0  0.0 

0.2  4 

X"  -  A*(l.  0  -  X**2)*X'  +  B*X  =  0.  0 


8.  0 


1.0 

1.0 

0.  0 

0.  1 

20.  0 

0.  01 

0.  0 

1 

0.  0 

0.  1 

20.  0 

1.  0 

0.  0 

1 

0.  0 

0.  1 

20.  0 

2.5 

3.  0 

1 

0.  0 

0.  1 

20.  0 

2.  0 

-3.  3 

1 

0.  0 

0.  1 

20.  0 

-2.  5 

3.  0 

1 

0.  0 

0.  1 

20.  0 

-2.  5 

-2.  0 
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Example   1:      X"  -  A*(l.  0  -  X**2)*X'  +  B*X  =  0.  0 


Fig.    5-1        X -SCALE  =  1.0  units /inch 
Y-SCALE  =1.0  units  /inch 
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Example  2: 
X"  +  A*X*X'  +  B*X  =  0.  0 
WHERE: 

A  =  1.0 

B  =  1.  0 

PLOT  PARAMETERS: 
XSCALE  =  1.0 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes     =     0.2 
Number  of  slopes  per  inch 


YSCALE 

= 

1.0 

YCENTER 

= 

0.  0 

YSIZE 

— 

8.  0 

=     4 


THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

X"  +  A#X*X'  +  B*X    =    0.0 

1.0  0.  0  6.  0  1.0 

0.2  4 

X"  +  A*X#X'  +  B*X  =  0.  0 


0.  0 


8.  0 


1.0 

1.  0 

0.  0 

0.  02 
1.0 

15.  0 

1 

0.  0 

0.  02 
2.  0 

15.  0 

1 

0.  0 

0.  02 
3.  0 

15.  0 

1 

0.  0 

0.  02 

15.  0 

2.5 

-2.  0 

1 

0.  0 

0.  02 

15.  0 

2.5 

-3.  5 
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Example  2:    X"  +  A*X*X'  +  B*X    =    0.  0 
A  =  1.  0,       B  =  1.  0 


N  N  \  \   I   \   \    \ 


Fig.      5-2        X-SCALE  =   1.  0  units /inch 
Y-SCALE  =  1.  0  units /inch 
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SCALE 

= 

1.0 

CENTER 

= 

0.  0 

SIZE 

= 

8.0 

Example  3: 

X"  +J1.  0  -  X**2)*X'  +  X  =  0.  0 

PLOT  PARAMETERS: 

XSCALE  =  1.0 

XCENTER        =  0.  0 

XSIZE  =  6.0 

Size  of  slopes     =     0.2 

Number  of  slopes  per  inch  =  4 

The  input  cards  used: 

NELSON,    H.    G. 

X"  +  (1.  0  -  X**2)*X'  +  X  =  0.  0 

1.0  0.0  6.0  1.0  0.0  8.0 

0.2  4 

X"  +  (1.  0  -  X**2)*X'  +  X  =  0.  0 

0.  0 

0.0  0.01  8.0 

1.0 

1 

0.0  0.01  10.0 

1.5 

1 

0.0  0.01  8.0 

2.0 

1 

0.  0  0.  01  8.  0 

2.25 

1 

0.  0  0.  01  8.  0 

-2.  05 
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Example  3:    XM  +  (1.  0  -  X**2)*X'  +  X  =  0.  0 


Fig.    5-3    X-SCALE 
Y-SCALE 


1.  0  units /inch 
1.  0  units /inch 
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Example  4: 
X"  +  X'**2  +  X  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  8.  0 

Size  of  slopes  —    0.2 
Number  of  slopes  per  inch  =  4 


YSCALE 

= 

1.0 

Y CENTER 

= 

0.  0 

YSIZE 

- 

6.0 

THE  INPUT   CARDS  USED: 


NELSON,    ! 

H.    G. 

X"  +  X' 

**2 

:  +  x  = 

0.  0 

2.  0 

0. 

0 

8.  0 

0.2 

4 

X'»  +  X' 

:  +  x  = 

0.  0 

0.  0 

0.  0 

0. 

05 

40. 

0 

-4.  0 

2. 

5 

1 

0.  0 

0. 

05 

40. 

0 

-7.  0 

1 

0.  0 

0. 

05 

40. 

0 

1.  0 

3. 

0 

1 

0.  0 

0. 

05 

40. 

0 

-5.  0 

1 

0.  0 

0. 

05 

40. 

0 

-3.  0 

1 

0.  0 

0. 

05 

40. 

0 

-1.0 

1 

0.  0 

0. 

05 

40. 

0 

0.4 

2. 

5 

1.  0 


0.  0 


6.  0 
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Example  5: 

X"  +  A*X  +  B*X**3  =  0.  0 
WHERE: 

A  =  -1.0 

B  =     0.25 

PLOT  PARAMETERS: 

XSCALE  =  1.0  YSCALE  =  1.0 

XCENTER       =  0.  0  YCENTER        =  0.  0 

XSIZE  =  8.  0  YSIZE  =  6.  0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

X"  +  A*X  +  B*X**3  =  0.  0 
1.0  0.0  8.0  1.0  0.0  6.0 

0.2  4 

X"  +  A*X  +  B*X**3  =  0.  0 

-1.0  0.25 

0.0  0.02  18.0 

1.0  1.2 

1 

0.0  0.02  18.0 

1.0  2.0 

1 

0.0  0.04  30.0 

0.  0  0.  03 

1 

0.0  0.02  18.0 

0.3 

1 

0.0  0.02  18.0 

-0.  3 

1 

0.0  0.02  18.0 

1.0 

1 

0.0  0.02  18.0 

-1.  0 

1 

0.0  0.02  18.0 

-1.  5 

0.0  0.02  18.0 

1.5 
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Example  6: 
X11  +  (1.  0  -  ABS(X))*X'  +  X  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  1.0 

XCENTER        =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

1.0 

Y CENTER 

= 

0.0 

YSIZE 

— 

8.  0 

THE  INPUT  CARDS  USED: 

NELSON,   H.    G. 

X"  +  (1.0-  ABS(X))*X*  +  X  =  0.  0 
1.0  0.  0  6.  0  1.0 

0.2  4 

X"  +(1.0-  ABS(X))*X'  +  X  =  0.  0 


15.  0 
15.  0 
15.  0 
15.  0 
15.  0 


0.  0 


8.  0 


0.  0 

0.  0 

0.  02 

1.0 

1 

0.  0 

0.  02 

1.5 

1 

0.  0 

0.  02 

2.  0 

1 

0.  0 

0.  02 

2.  8 

1 

0.  0 

0.  02 

-2.5 
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Example  6:    X"     +(1.0-  ABS(X))*X'  +  X  =  0.  0 


Fig.    5-6      X-SCALE  =  1.  0  units/inch 
Y-SCALE  =  1.  0  units /inch 
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YSCALE 

= 

2.  0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

8.  0 

Example  7: 
X"  +  X'    +  X*ABS(X)  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.  0 

XCENTER        =  0.0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

X"  +  X'  +  X*ABS(X)  =  0.  0 
2.0  0.0  6.0  2.0  0.0  8.0 

0.2  4 

X"  +  X1  +  X*ABS(X)  =  0.  0 

0.  0 

0.0  0.01  8.0 

3.8 

1 

0.0  0.01  8.0 

-3.  0 

1 

0.  0  0.  01  8.  0 

5.0 

1 

0.0  0.01  8.0 

-5.  5 
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Example  7:    X"  +  X'  +  X*ABS(X)  =  0.  0 


Fig.    5-7      X-SCALE  =  2.  0  units/inch 
Y-SCALE  =  2.0  units/inch 
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C.      VARIOUS  TYPES  OF  PENDULUM  EXAMPLES 

Example  8: 

Linear  approximation  to  the  pendulum,    g/1  =  1 
No  Damping 

X"  +  X  =  0 

PLOT  PARAMETERS: 

XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

2.  0 

Y CENTER 

= 

0.  0 

YSIZE 

= 

8.  0 

THE  INPUT  CARDS  USED: 
NELSON,    H.    G.         PENDULUM 

X"  +  X  =  0 
2.0  0.  0  6.  0  2.0  0.  0 

0.2  4 

X"  +  X  =  0 


8.  0 


0.  0 
2.  0 
1 

0.  0 
1.0 
1 
0.  0 

4.  0 
1 

0.  0 

5.  5 


0.  02 


0.  02 


0.  02 


0.  02 


10.  0 
10.  0 
10.  0 
10.  0 
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Example  8:    X"  +  X  =  0 
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Fig.    5-8      X-SCALE  =  2.  0  units/inch 
Y-SCALE  =  2.0  units/inch 
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Example  9: 
Pendulum  without  damping,    g/1  =  1 
X"  +  SIN(X)  =  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  8.  0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 
THE  INPUT  CARDS  USED: 

NELSON,    H.    G.  PENDULUM 

X"  +  SIN(X)  =  0 
2.0  0.0  8.0  2.0  0.0 

0.  2  4 

X"  +  SIN(X)  =  0 


YSCALE 

~ 

2.  0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

6.  0 

6.  0 


0.  0 

0.  02 

10.  0 

-6.0 

1.  o 

1 

0.  0 

0.  02 

15.  0 

-6.  0 

2.  0 

1 

0.  0 

0.  02 

10.  0 

-6.  0 

3.  0 

1 

0.  0 

0.  02 

10.  0 

-6.  0 

4.  0 

1 

0.  0 

0.  02 

10.  0 

0.  0 

1.  0 

1 

0.  0 

0.  02 

10.  0 

6.0 

1.  o 

1 

0.  0 

0.  02 

15.  0 

6.  0 

-2.  0 

1 

0.  0 

0.  02 

10.  0 

6.  0 

-3.  0 

1 

. 

0.  0 

0.  02 

10.  0 

6.0 

-4.  0 

68 


69 


Example  10: 
Pendulum  with  velocity  damping. 
X"  +  A*X'+  SIN(X)  =  0 
WHERE: 
A  =  0.  3 

PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER       =  0.  0 

XSIZE  =  8.  0 

Size  of  slopes  =  0.  1 
Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

2.0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

6.0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G.  PENDULUM 

X"  +  A*X'  +  SIN(X)  =  0 
2.  0  0.  0  8.  0  2.  0  0.  0 

0.  1  4 

X"  +  A*X'  +  SIN(X)  =  0 


6.  0 


0.  3 

0.  0 

0.  05 

25.  0 

-7.0 

3.  0 

1 

0.  0 

0.  05 

25.  0 

-7.  0 

5.  0 

1 

0.  0 

0.  05 

25.  0 

-7.  0 

2.  0 

1 

0.  0 

0.  05 

25.  0 

6.  0 

-5.  0 
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Example  1 1 : 
Pendulum  with  velocity  squared  damping 
X"  +  A*X'*ABS(X')  +  SIN(X)  =  0 
WHERE: 
A  =  1.  0 

PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  8.  0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

2.  0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

6.0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G.  PENDULUM 

X"  +  A*X'*ABS(X')  +  SIN(X)  =  0 

2.  0  0.  0  8.  0  2.  0 

0.2  4 

X"  +  A*X!*ABS(X')  +  SIN(X)  =  0 


20.  0 


1.0 

0.  0 

0.  03 

1.0 

4.  0 

1 

0.  0 

0.  03 

-4.  0 

4.  0 

1 

0.  0 

0.  03 

-7.  0 

4.  0 

1 

0.  0 

0.  03 

5.  0 

-4.  0 

1 

0.  0 

0.  03 

2.0 

4.  0 

1 

0.  0 

0.  03 

-2.  5 

-4.  0 

20.  0 


20.  0 


20.  0 


20.  0 


20.  0 


0.  0 


6.  0 
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D.      SATURATED  SYSTEMS 

Example  12: 

A  saturated  servo  system.     Saturation  occurs  when  absolute  error 

signal  exceeds   0.  2. 

IF(X.  LT.  -0.  2)    THEN  X"  +  0.  2*X'  +  X  =  0.  2 
IF(X.  GT.  0.  2)      THEN  X"  +  0.  2*X'  +  X  =  -0.  2 
X"  +  0.2*X'  +  X  =  0.  0 

PLOT  PARAMETERS: 

XSCALE  =  0.3  YSCALE  =  0.2 

XCENTER        =  0.  0  YCENTER        =  0.  0 

XSIZE  =  6.  0  YSIZE  =  8.  0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 


THE  INPUT  CARDS  USED: 

NELSON,    H.    G.      0902 

SATURATED  SERVO  SYSTEM 
0.  3  0.  0  6.0  0.  2  0.  0  8.0 

0.  2  4 

IF(X.  LT.  (1.  0-2.  0)*0.  2)    THEN  X"  +  0.  2*X'  +  X 

IF(X.  GT.  0.  2)    THEN  XM  +  0.2*X'  +  X  =  -0.2 

X"  +  0.  2*X'  +  X  =  0.  0 


=  0.2 


0.  0 
0.7 


0.  05 


30.  0 
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Example  12:  Saturated  Servo  System 

/  •/  /  \  / 

/  /  /  v  / 

/  /  / 


'  /  /  ih  i  i  : 


/  r  y  /  I  I 


\\\  \\\\ 

Fig.    5-12      X-SCALE  =  0.  3  units  /inch 
Y-SCALE  =0.2  units /inch 


■      /    /     I 

/  /  >  /  /  i 


75 


E.      RELAY  SYSTEMS 

Example  13: 

IDEAL  RELAY 

X"  +  0.  3*X'  +  1.  5*SIGN(X)  =  0.  0 

PLOT  PARAMETERS: 

XSCALE  =  1.0  YSCALE 

XCENTER        =  0.  0  YCENTER 

XSIZE  =  6.0  YSIZE 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 


0.6 
0.  0 

8.  0 


THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 
IDEAL  RELAY 
1.0  0.0  6.0  0.6 

0.2  4 

X"  +  0.  3*X'  +  1.  5*SIGN(X)  =  0.  0 


0.  0 


8.  0 


0.  0 
2.8 


0.  1 


85.  0 
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Example  13:      X"  +  0.  3*X'  +  1.  5*SIGN(X)  =  0.  0 


\    \ 
\ 
\    \ 


\    \    \    \    \ 


Fig.    5-13      X-SCALE  =  1.  0  units  /inch 
Y-SCALE  =   0.  6  units /inch 
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Example  14: 
Ideal  Relay  with  Rotated  Switching  Line 
X"  +  0.  2*X'  +  1.  0*SIGN(4.  0*X  -  X!)  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.  0 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

= 

1.0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

8.0 

THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 

IDEAL  RELAY  WITH  ROTATED  SWITCHING  LINE 
2.0  0.0  6.0  1.0  0.0  8.0 

0.2  4 

X"  +  0.  2*X'  +  1.  0*SIGN(4.  0*X  -  X')  =  0.  0 


0.  0 
0.  0 

5.  5 


0.  1 


80.  0 
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Example  14:    X"  +  0.  2*X'  +  1.  0*SIGN(4.  0*X  -  X1)  =  0.  0 


Fig.    5-14      X-SCALE  =  2.  0  units /inch 
Y-SCALE  =   1.  0  units /inch 
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Example  15: 
Ideal  Relay  with  Rotated  Switching  Line 
X"  +  0.  2*X'  +  1.  0*SIGN(4.  0*X  +  X1)  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER       =  0.  0 

XSIZE  =  6.  0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

0.8 

Y CENTER 

= 

0.  0 

YSIZE 

~ 

8.0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

IDEAL  RELAY  WITH  ROTATED  SWITCHING  LINE 
2.0  0.0  6.0  0.8  0.0  8.0 

0.2  4 

X"  +  0.  2*X'  +  1.  0*SIGN(4.  0*X  +  X1)  =  0.  0 


0.  0 
0.  0 
5.  5 


0.  1 


20.  0 


80 


Example  15:      X"  +  0.2*X'  +  1.  0*SIGN(4.  0*X  +  X1)  =  0.  0 


\    \    \    >    \ 

\  \  I 
\\\U 

Fig.    5-15      X-SCALE  =  2.  0  units/inch 
Y-SCALE  =   0.  8  units /inch 
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Example  16: 
Ideal  Relay  with  Rotated  Switching  Line 
X"  +  0.2*X'  +  1.  0*SIGN(2.  0*X  +  X')  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

2 

1.  0 

Y CENTER 

= 

0.  0 

YSIZE 

— 

8.  0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

IDEAL  RELAY  WITH  ROTATED  SWITCHING  LINE 
2.0  0.0  6.0  1.0  0.0  8.0 

0.2  4 

X"  +  0.  2*X'  +  1.  0*SIGN(2.  0*X  +  X')  =  0.  0 


0.  0 
0.  0 
5.  5 


0.  1 


20.  0 
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Example  16:    X"  +  0.  2*X'  +  1.  0*SIGN(2.  0*X  +  X1)  =  0.  0 


Fig.    5-16      X-SCALE  =  2.  0  units  /inch 
Y-SCALE  =   1.  0  units /inch 
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Example  17: 
Ideal.  Relay  with  Rotated  Switching  Line 
XM  +  0.2*X'  +  1.  0*SIGN(2.  0*X  -  X')  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.  0 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

— 

0.  8 

Y CENTER 

= 

0.  0 

YSIZE 

- 

8.  0 

THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 

IDEAL  RELAY  WITH  ROTATED  SWITCHING  LINE 
2.0  0.0  6.0  0.8  0.0  8.0 

0.2  4 

X"  +  0.  2*X'  +  1.  0*SIGN(2.  0*X  -  X')  =  0.  0 


0.  0 
0.  0 

5.  5 


0.  1 


80.  0 
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Example  17: 


X"  +  0.  2*X*  +  1.  0*SIGN(2.  0*X  -  X')  =  0.  0 

I 


\    \    \ 

\     \      \       \      \      \       \       * 

Fig.    5-17      X-SCALE  =  2.  0  units/inch 
Y-SCALE  =  0.  8  units /inch 
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F.      DEAD  ZONE 

Example   18: 

IF(X.  LT.  -A.  OR.  X.  GT.  A)  THEN  X"  +  B*X'  +  C*X  =  C*A*SIGN(X) 
X"  +  B*X'  =  0.  0 

WHERE: 

A  =  0.  3 
B  =  0.  2 
C  =  1.0 

PLOT  PARAMETERS: 

XSCALE  =  0.  5 

XCENTER       =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 


YSCALE 

~ 

0.3 

Y CENTER 

= 

0.  0 

YSIZE 

- 

'    8.  0 

THE  INPUT   CARDS  USED:  ; 


NELSON,    H.    G. 

DEAD  ZONE 
0.  5  0.  0  6.  0  0.  3  0.  0  8.  0 

0.  2  4 

IF(X.  LT.  -A.  OR.X.  GT.A)  THEN  X"  +  B*X'  +  C*X 

X"  +  B*X'  =   0.  0 


C*A*SIGN(X) 


0.3 

0.2 

1.  0 

0.  0 

0.  1 

50.  0 

0.  0 

1.3 

1 

0.  0 

0.  1 

50.  0 

1.  0 
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Example  18:      Dead  Zone 

IF(X.  LT.  -A.  OR.X.  GT.A)  THEN  X"  +  B*X'  +  C*X  =  C*A*SIGN(X) 

X"  +  B*X'  =  0.  0 

A-=  0.  3,    B  =  0.2,    C  =  1.0 


X  /  /  / 
/  X  /  /  / 


Fig.    5-18      X-SCALE  =  0.  5  units /inch 
Y-SCALE  =  0.  3  units /inch 
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Example   19:     Dead  Zone 

IF(ABS(X).  LT.  2.  0)  THEN  X"  +  0.  2*X'  =  0.  0 
X"  +-0.  2*X'  +  X  =  0.  0 


PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER       =  0.  0 

XSIZE  =  8.  0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

— 

.     2.  0 

YCENTER 

= 

0.  0 

YSIZE 

= 

6.  0 

THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 
DEAD  ZONE 

2.0  0.0  8.0  2.0  0.0  6.0 

0.2  4 

IF(ABS(X).  LT.  2.  0)  THEN  X"  +  0.  2*X'  =  0.  0 

X"  +  0.  2*X'  +  X  =  0.  0 


0.  0 
0.  0 
6.  0 


0.  05 


35.  0 
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G.      A  SECOND  ORDER  LINEAR  SYSTEM 

Example  2  0: 

The  following  three  plots  demonstrating  the  effect  of  various 

damping  were  the   result  of  the  input  cards  listed  below. 

X"  +  2.  0*A*B*X'  +  B**2*X  =  0.  0 

"Where: 

A  =  is  the  natural  frequency 
B  =  is  the  damping  factor 

PLOT  PARAMETERS: 

XSCALE  =  2.0  YSCALE  =  2.0 

XCENTER        =  0.  0  YCENTER        =  0.  0 

XSIZE  =  8.0  YSIZE  =  6.0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 


THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 

X"  +2.  0*A*B*X'  +  B**2*X  =  0.  0 
2.0  0.  0  8.  0  2.  0  0.  0 

0.2  4 

X"  +  2.  0*A*B*X'  +  8**2 *X  =  0.  0 


6.  0 


1.0 

1.  0 

0.  0 

0.  1 

7.  0 

1 

0.  0 

0.  1 

-7.  0 

1 

0.  0 

0.  1 

3.  5 

1 

0.  0 

0.  1 

-3.  5 

2 

0.  5 

1.  o 

0.  0 

0.  1 

7.  0 

1 

10.  0 
10.  0 
15.  0 
10.  0 

15.  0 
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0.  0 
-7.  0 
1 

0.  0 
3.5 
1 

0.  0 
-3.  5 
2 

0.  1 
0.  0 
5.0 


0.  1 


0.  1 


0.  1 


1.0 
0.  1 


15.  0 


15.  0 


15.  0 


80.  0 
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H.      A  FIFTH  ORDER  SYSTEM 

Example  21: 

IF(T     .GT.     0.  0015)  SKIP  3 

VI  =  3000511. 6  +  A 

V2  =  285120.  0  +  A*3511.  0 

V3  =  A*38500.  0 

X,IMI  +  122.  2*X,m  +  50835.  8*X,M  +  V1*XM  +  V2*X'  +  V3*X    =  0.  0 


WHERE: 

G  =  30212. 0 

PLOT  PARAMETERS: 

XSCALE 

2.  0 

YSCALE 

20.  0 

XCENTER 

0.  0 

Y CENTER 

0.  0 

XSIZE 

6.0 

YSIZE 

8.  0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 
FIFTH  ORDER  SYSTEM 
2.0  0.  0  6.  0  20.  0  0.  0  8.  0 

IF(T.GT.    0.  0015)  SKIP  3 

VI  =  3000511.  6  +  A 

V2  =  285120.  0  +  A*3511.  0 

V3  =  A*38500.  0 

X,MM  +  122.  2*XMn  +  50835.  8*X"'  +  V*X"  +  V2*X'     +   V3*X  =  0.  0 


0.4 


30212. 

0 

0.  0 

0. 

001 

-3.  0 

1 

0.  0 

0. 

001 

2.5 

1 

0.  0 

0. 

001 

-5.  1 

1 

0.  0 

0. 

001 

4.6 

0.4 


0.4 


0.4 
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Example  21:      A  FIFTH  ORDER  SYSTEM 

IF(T    .GT.     0.0015)    SKIP  3 

VI  =  3000511. 6  +  A 

V2  =..285120.  0  +  A*3511.  0 

V3  =  A*38500.  0 

X"1"  +  122.  2*X"M  +  50835.  8*X'"  +  V1*X"  +  V2*X'  +  V3*X 


=  0.  0 


Fig.    5-21      X-SCALE  =  2.  0  units /inch 
Y-SCALE  =  20.  0  units/inch 
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I.       HYSTERESIS  AND  DEAD  ZONE 

Example  22: 

RELAY  WITH  DEAD  ZONE  AND  HYSTERESIS 

THE  SYSTEM  DESCRIPTION  EQUATIONS  AND/OR  PROGRAM: 

V3  =  T 

IF(ABS(X)    .  LE.     A)V1  =  0.  0 

IF(ABS(X)    .  GT.     B)V1  =  SIGN(X) 

X"  +  0.2*X'  +  VI  =  0.  0 

WHERE: 

A  =  0.  8 

B  =   1.2 

The  first  statement  is  used  to  suppress  the  slopes. 
PLOT  PARAMETERS: 


XSCALE 

= 

1.0 

YSCALE 

= 

0.  5 

XCENTER 

= 

0.  0 

YCENTER 

= 

0.  0 

XSIZE 

— 

6.0 

YSIZE 

- 

8.  0 

THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 

RELAY  WITH  DEAD  ZONE  AND  HYSTERISES 
1.0  0.0  6.0  0.5  0.0  8.0 

V3  =  T 

IF(ABS(X).  LE.A)  VI  -  0.  0 
IF(ABS(X).  GT.  B)  VI  =  SIGN(X) 
X"  +  0.  2*X'  +  VI  =  0.  0 


30.  0 


0.  8 

1.2 

0.  0 

0.  04 

-2.  8 
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Example  22:      RELAY  WITH  DEAD  ZONE  AND  HYSTERESIS 


Fig.    5-22      X-SCALE  =  1.  0  units/inch 
Y-SCALE  -   0.  5  units /inch 
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J.       NON-LINEAR  RESTORING  FORCES 

Example  23: 

SECOND  ORDER  SYSTEM  WITH  A  RESTORING  FORCE  CORRESPONDING 
TO  A  NON- LINEAR  SPRING 

THE  SYSTEM  DESCRIPTION  EQUATION: 

X"  +  0.  9*X'  +  (X  +  A*X**3)  =  0.  0 

WHERE: 

A  =  0.  2 

PLOT   PARAMETERS: 

XSCALE  =  2.  0 

XCENTER        =  0.  0 

XSIZE  =  6.0 

Size  of  slopes  =  0.  2 

Number  of  slopes  per  inch  =  4 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

SECOND  ORDER  SYSTEM  WITH  NON-LINEAR  SPRING 
2.0  0.0  6.0  2.0  0.0  8.0 

0.2  4 

X"  +  0.  9*X'  +  (X  +  A*X**3)  =  0.  0 


YSCALE 

— 

2.  0 

YCENTER 

= 

0.  0 

YSIZE 

— 

8.  0 

0.2 

0.  0 

0.  02 

15.  0 

-4.  0 

2 

-0.  1 

0.  0 

0.  02 

15.  0 

-3.  0 

1 

0.  0 

0.  02 

15.  0 

1.0 

-7.  0 

1 

0.  0 

0.  02 

15.  0 

-5.  5 

5.  0 

1 

0.  0 

0.  02 

15.  0 

-5.2 

5.  0 

1 

0.  0 

0.  02 

15.  0 

5.  0 

-5.  0 
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Example  23(a):      X"  +  0.  9*X*  +  (X  +  A*X**3)  =  0. 
A  =  -0.  1  (gives  a  "soft"  spring) 

\    \    \  \  \  \  \  ^  X  X  X  X  X  \  \  \ 


\    \    \  \  X  \  X  X  X  X  X 


v 


Fig.    5-23(a) 


X -SCALE 
Y-SCALE 


=  2.0  units  /inch 
=  2.0  units  /inch 


N  \  \  % 

••.  \ 

\ 

wx1 

N    \ 

\ 

\  \  \  s' 

h    \ 

\ 
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Example  23(b):      X"  +  0.  9*X'  +  (X  +  A*X**3 )  =  0.  0 


a  =  o. ; 

l  (gives  a  "1 

iard' 

1  spring) 
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Fig.    5-23(b)         X-SCALE  =  2.0  units/inch 
Y-SCALE  =  2.0  units/inch 
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YSCALE 

= 

1. 

0 

YCENTER 

= 

0. 

0 

YSIZE 

s 

8. 

0 

Example  24: 

A  FORCED  SYSTEM  WITH  A  NON- LINEAR  RESTORING  FORCE 

THE  SYSTEM  DESCRIPTION: 

X"  +  0.  9*X!  +  (X  +  0.  2*X**3)  =  A*COS(B*T) 

WHERE: 

A  =  3.  0 
B  =  0.  7 

PLOT  PARAMETERS: 

XSCALE         =         1.0 

XCENTER    =        0.0 

XSIZE  =        6.  0 

THE  INPUT   CARDS  USED: 

NELSON,    H.    G. 

FORCED  NON- LINEAR  SECOND  ORDER  SYSTEM 
1.0  0.  0  6.  0  1.0  0.  0  8.  0 

X"  +0.  9*X'  +  (X  +  0.  2#X**3)  =  A*COS(B*T) 

3.0  0.  7 

0.0  0.04  30.0 

0.  0  0.  0 

3 

NELSON,    H.    G. 

FORCED  NON- LINEAR  SECOND  ORDER  SYSTEM 
2.  0  0.  0  6.  0  2.0  0.  0  8.  0 

X"  +  (X  +  0.  2*X**3)  =  A*COS(B*T) 

5.0  0.  8 

0.  0  0.  04  30.  0 

0.  0  0.  0 

2 

5.  0  1.0 

0.  0  0.  04  30.  0 

0.  0  0.0 
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Example  24(a):      X"  +0.  9*X'  +  (X  +  0.  2*X**3)  =  A*COS(B*T) 
A  =  3.  0,       B  =  0.  7 
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1  "  *1     J 
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a 

Fig.    5-24(a)      X-SCALE  =  1.  0  units /inch 
Y-SCALE  =  1.  0  units /inch 
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Example  24(b):      XM  +  0.  0*X'  +  (X  +  0.  2*X**3)  =  A*COS(B*T) 
A=  5.  0,     B  =  0.  8 


Fig.    5-24(b)      X-SCALE 
Y-SCALE 


1.  0  units  /inch 
1.  0  units  /inch 
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Example  24(c):      X"  +  0.  0*X'  +  (X  +  0.  2*X**3)  =  A*COS(B*T) 
A  =  5.  0,     B  =  1.  0 


Fig.    5-24(c)      X-SCALE  =   1.  0  units /inch 
Y-SCALE  =   1.  0  units /inch 
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K.      BANG-BANG  SYSTEMS 

Example  25. 
A  SECOND  ORDER  BANG-BANG  SYSTEM 
SYSTEM  EQUATION: 
X"  +  SIGN(X  +  0.  5*X'*ABS(X'))  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  8.0 

Size  of  slopes  =  0.  2 
Number  of  slopes  per  inch  =  4 


YSCALE 

= 

1.0 

Y  CENTER 

= 

0.  0 

YSIZE 

- 

6.  0 

THE  INPUT  CARDS  USED: 

NELSON,    H.    G. 

SECOND  ORDER  BANG-BANG 
2.0  0.0  8.0  1.0  0.0 

0.2  4 

X"  +  SIGN(X  +  0.  5*X,*ABS(X1))  =  0.  0 


6.  0 


0.0  0.02  15.0 

-7.  0 

1 

0.0  0.02  15.0 

7.  0 
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Example  26. 
A  THIRD  ORDER  BANG-BANG  SYSTEM 
SYSTEM  EQUATIONS: 
VI  =  SIGN(X'  +  0.  5*X"*ABS(X")) 
V2  =  X  +  (1.  /3.  )*X"**3  +  Vl*Xn#X' 
V3  =  VI  *(0.  5*X"**2  +  V1*X')**1.  5 
X"1  +  SIGN(V2  +  V3)  =  0.  0 
PLOT  PARAMETERS: 
XSCALE  =  2.0 

XCENTER        =  0.  0 

XSIZE  =  6.  0 


The  input  cards  used: 
NELSON,    H.    G. 
THIRD  ORDER  BANG-BANG  SYSTEM 

2.0  0.0  6.0  1.0  0.0  8.0 

VI  =  SIGN(X'  +  0.  5*X"*ABS(X")) 
V2  =  X  +  (1.  /3.  )*X"**3  +  V1*XM*X 
V3  =  V1*(0.  5#XM**2  +  V1*X')**1.  5 

X'"  +  SIGN(V2  +  V3)  =  0.  0 


YSCALE 

— 

1.0 

Y CENTER 

= 

0.  0 

YSIZE 

~ 

8.  0 
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Example  26:      A  THIRD  ORDER  BANG-BANG  SYSTEM 
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Fig.    5-26      X-SCALE  =  2.  0  units /inch 
Y-SCALE   =   1.  0  units /inch 
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VI.    "CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

This  thesis  has  introduced  an  easy  to  use,    flexible,    and  powerful 
analysis  tool  for  the  study  of  non-linear  ordinary  differential  equations. 
The  primary  application  emphasis  has  been  in  relation  to  control 
system  problems  for  which  many  examples  have  been  presented.      It 
has  been  shown  that  the  use  of  an  automatic  built-in  compiler  coupled 
with  interactive  graphics   result  in  a  very  useful  analysis  program. 
For  a  large  number  of  quite  diverse  examples  it  has  been  demonstrated 
in  this  thesis  that  one  can  very  adequately  specify  a  problem  using 
data  cards  or  by  typing  on  an  interactive  graphics  screen. 

This  study  has  suggested  several  basic  ideas  for  further  research. 

1.  An  extension  of  the  program  following  the  development  presented 
in  Chapter  II  would  give  the  program  the  capability  of  solving  systems 
described  by  several  simultaneous  non-linear  differential  equations. 

2.  If  one  were  concerned  with  execution  efficiency,    it  would  be 
very  worthwhile  to  rewrite  the  interpreter  in  assembly  language.     The 
interpreter  is  called  four  times  for  each  solution  point  and  once  for 
each  slope  line.     Thus  the  slope  generation  and  solution  generation 
times  could  be  speeded  up  considerably  with  a  faster  interpreter.     The 
other  parts  of  the  program  are  used  relatively  infrequently  and  there- 
fore efficiency  is  not  an  important  factor. 

3.  It  would  be  worthwhile  to  provide  a  choice  of  integration 
schemes  to  the  user  in  addition  to  the  present  fourth  order  Runge- 
Kutta . 
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4.  A  built-in  compiler  could  be  used  in  several  other  applications. 
One  of  the  very  obvious  areas  is  that  of  optimization.      For  example, 

a  very  useful  directed  search  algorithm  called  "DIRECT"  is  available 
on  the  XDS  9300.     The  function  to  be  minimized  is  incorporated  into  a 
function  subroutine  which  must  be  supplied  by  the  user.     A  modification 
of  the  compiler  and  interpreter  presented  in  this  thesis  could  be  used 
as  a  very  powerful  interface  for  describing  the  function  and  its 
boundaries. 

5.  It  would  be  very  worthwhile  developing  methods  to  assist  in 
the  analysis  of  systems  higher  than  second  order.      In  general,    the 
presentation  of  the  grid  of  slopes  is  not  worthwhile  for  systems  of 
higher  than  second  order  since  slope  lines  are  in  a  space  which  has  a 
dimension  equal  to  the  order  of  the  differential  equation. 
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APPENDIX    A 
REQUIRED  JOB  CONTROL  LANGUAGE  FOR  BATCH  VERSION 


1.  When  using  the  compiled  binary  object  deck  use  the  following 
procedure: 

//      JOB  CARD  -Green- 

//      EXEC  FORTLG,  REGION.  GO=160K 

//LINK.  SYSIN  DD  * 

BINARY  OBJECT  DECK 

//GO.  FT06F001     DD    S  YSOUT=A,  SPACE  =  (CYL,  (6)) 
//GO.SYSPLOT  DD  SYSOUT  =  C,  SPACE=(TRK,  (1,  10)) 
//GO.SYSPLOTS  DD  UNIT=S  YSDA,  SPACE=(TRK,  (2,  6)) 
//GO.  SYSIN  DD  * 

DATA 

/*  -Orange- 

2.  If  using  the  source  deck  use  the  following  procedure: 

//        JOB     CARD  -Green- 

//        EXEC  FORTCLGP,  REGION.  FORT  =  150K,  REGION.  GO=160K 
//FORT.  SYS  PRINT  DD  SYSOUT=A,  SPACE  =  (CYL,  (1,  1)) 
//FORT.  SYSIN  DD  * 

FORTRAN  SOURCE  PROGRAM 

/* 

//GO.  FT06F001  DD  SYSOUT=A,  SPACE(CYL,  (6)) 
//GO.SYSPLOT  DD  SYSOUT  =  C,  SPACE=(TRK,  (1,  10)) 
//GO.SYSPLOTS  DD  UNIT=S YSDA,  SPACE=(TRK,  (2,  6)) 
//GO.  SYSIN  DD  * 

DATA 

/*  -  Orange- 
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APPENDIX    B 

1.       Required  control  cards  for  Interactive  Graphics  Version. 

BOOT    CARD 

PATCH  DECK 

AAGT 
4  JOB 

L  ASSIGN  BI  =MTOA 

A  ASSIGN  GO=MTOA 

A  LOAD  XM,    MAP 

A  SEG  MAIN-(SINIT,    SSETUP2,  SCOMP,  (SGRID-SINTER),  (SSOLVE- 

ADATA  SINTER)) 
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PROGRAM  A 
LISTING  OF  BATCH  VERSION 
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